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ABSTRACT 



Energy-momentum is an important conserved quantity whose definition has 
been a focus of many investigations in general relativity. Unfortunately, there 
is still no generally accepted definition of energy and momentum in general 
relativity. Attempts aimed at finding a quantity for describing distribution of 
energy-momentum due to matter, non-gravitational and gravitational fields 
only resulted in various energy-momentum complexes (these are nontensorial 
under general coordinate transformations) whose physical meaning have been 
questioned. The problems associated with energy-momentum complexes re- 
sulted in some researchers even abandoning the concept of energy-momentum 
localization in favor of the alternative concept of quasi-localization. However, 
quasi-local masses have their inadequacies, while the remarkable work of 
Virbhadra and some others, and recent results of Cooperstock and Chang et 
al. have revived an interest in various energy-momentum complexes. Hence 
in this work we use energy-momentum complexes to obtain the energy dis- 
tributions in various space-times. 

We elaborate on the problem of energy localization in general relativity 
and use energy-momentum prescriptions of Einstein, Landau and Lifshitz, 
Papapetrou, Weinberg, and M0ller to investigate energy distributions in var- 
ious space-times. It is shown that several of these energy-momentum com- 
plexes give the same and acceptable results for a given space-time. This shows 
the importance of these energy-momentum complexes. Our results agree with 
Virbhadra's conclusion that the Einstein's energy-momentum complex is still 
the best tool for obtaining energy distribution in a given space-time. The 
Cooperstock hypothesis (that energy and momentum in a curved space-time 
are confined to the the regions of non- vanishing energy-momentum of matter 
and the non-gravitational field) is also supported. 
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Chapter 1 
Introduction 



The notions of energy-momentum together with conservation laws play a 
fundamental role in any physical theory. The importance of conservation of 

energy and momentum concepts was first clearly identified more than two 
centuries ago within the Newtonian mechanics of a closed system of mass 
points. In this case both the linear momentum and the sum of kinetic and 
potential energies (for conservative forces) are conserved. The connection be- 
tween these conservation laws and their invariance under Galilei group was 
only detected around the year 1900 (see in Havas[30]). In continuum mechan- 
ics of elastic bodies the formulation of energy conservation laws necessitated 
the introduction of the concept of strain energy. Maxwell and Poynting had 
to introduce the concepts of electric and magnetic energy densities in classi- 
cal electrodynamics to retain energy conservation. Einstein's famous result 
of the special theory of relativity (SR) that mass is equivalent to energy 
is a consequence of the requirement that the law of conservation of energy 
and momentum should hold in all inertial frames. As a matter of fact, this 
process of introducing new forms of energy in order to retain conservation 
laws characterizes the whole development of physics. However, this practice 
involving the introduction of new kinds of energy ran into serious difficulties 
with the arrival of the general theory of relativity (GR). The main difficulty 
is with the expression defining the gravitational field energy part. 

Now, after more than eighty years of the success story of the theory of 
General Relativity (this includes, amongst other things, verification of the 
deflection of light by the Sun, the perihelion advance of Mercury, the gravita- 
tional red shift of light, and discoveries of quasars, cosmic fireball radiation, 
pulsars. X-ray sources that might contain black holes, and the present inter- 



1 



est in the imminent detection of gravitational waves), there is still no gen- 
eral agreement on the definition of energy, and more generally, of conserved 
quantities associated with the gravitational field. This dilemma in GR is 
highlighted in an important paper by Penrose[53] in the following way: "It is 
perhaps ironic that energy conservation . . . which now has found expression 
in the (covariant) equation 

VaT'^' = 0, (1.1) 

. . . should nevertheless have found no universally applicable formulation, within 
Einstein's theory, incorporating the energy of gravity itself". Indeed, Ein- 
stein's search for his generally covariant field equations was not only guided 
by the principle of equivalence but also by conservation laws of energy- 
momentum. Conservation laws of energy-momentum played a major role in 
the development of Einstein and Grossmann's so-called Entwurf theory [49]. 
Although the field equations of the Entwurf theory were only of limited co- 
variance, the theory had all the essential features of Einstein's final theory 
of GR. Einstein and Grossmann had considered the use of the Ricci tensor 
in deriving almost covariant field equations but had rejected these equations 
because of misconceptions which were based on Einstein's earlier work on 
static gravitational fields (for details see in Norton [49]). 

In his[23] derivation of the generally covariant gravitational field equa- 
tions, Einstein formulated the energy-momentum conservation law in the 
form: 

^(x/=^(7;>+V)) = o- (1-2) 

With Tj^ representing the stress energy density of matter^, Einstein iden- 
tified tj^ as representing the stress energy density of gravitation. Einstein 
noted that t^* was not a tensor, but concluded that the above equations 
(1.2) hold good in all coordinate systems since they were directly obtained 
from the principle of general relativity. The choice of a nontensorial quan- 
tity to describe the gravitational field energy immediately attracted some 
criticism. Levi-Givita not only attacked Einstein's use of a pseudotensor 
quantity (which is only covariant under linear transformations) in describing 
the gravitational field energy, but also suggested an alternative gravitational 
energy tensor which required that Einstein's field equations be also inter- 
preted as conservation laws (for details see in Gattani and De Maria [11], 

-^"Matter" includes the energy contribution of all non-gravitational fields. 
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and in Pauli [52] ). Both Schrodinger and Bauer came up with counter ex- 
amples to Einstein's choice of a nontensor. Schrodinger showed that by a 
suitable choice of a coordinate system the pseudotensor of a Schwarzschild 
solution vanishes everywhere outside the Schwarzschild radius. Bauer's ex- 
ample illustrated that a mere introduction of polar coordinates instead of 
quasi- Cartesian coordinates into a system of inertia without matter present 
would create a nonvanishing energy density in space. By resorting to the 
equivalence principle and physical arguments, Einstein vigorously defended 
the use of his pseudotensor to represent gravitational field[ll]. The prob- 
lems associated with Einstein's pseudotensor resulted in many alternative 
definitions of energy, momentum and angular momentum being proposed 
for a general relativistic system (see Aguirregabiria et al [1] and references 
therein) . 

The lack of a generally accepted definition of energy distribution in curved 
space-times has led to doubts regarding the idea of energy localization. The 
ambiguity in the localization of energy is not a new physics problem, pe- 
culiar to the theory of GR, but is also present in classical electrodynamics 
(Feynmann et al[25]. In GR there is a dispute with the importance of non- 
tensorial energy-momentum complexes^ whose physical interpretation has 
been questioned by a number of scientists, including Weyl, Pauli and Ed- 
dington (see in Chandrasekhar and Ferrari [15]). There are suspicions that, 
in a given space-time, different energy distributions would be obtained from 
different energy-momentum complexes. However, Virbhadra and co-workers 
investigated several examples of particular space-times (the Kerr-Newman, 
the Einstein-Rosen, and the Bonnor-Vaidya) and found that different energy- 
momentum complexes give the same energy distribution for a given space- 
time. Several energy-momentum complexes have been shown to coincide for 
any Kerr-Schild class metric [1]. In this thesis we are extending the work of 
Virbhadra and co-workers by considering further space-times in showing that 
the energy-momentum complexes are useful in obtaining meaningful energy 
distribution in a given geometry. In the rest of the chapter we give a brief 
review of both the theory of special relativity and general relativity so as to 
establish both the notation and terminology which will be used in the rest of 
our discussion. We also give a list of formulas needed in our later discussion. 

^We use the term energy-momentum complex for one which satisfies the local conserva- 
tion laws and gives the contribution from the matter (including all nongravitational fields) 
as well as the gravitational field. 
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1.1 Tensor equations 



Einstein's 1905 paper on the special theory of relativity which demanded a 
complete change of attitude towards space and time prompted Minkowski 
to declare with great elegance, in his famous 1908 speech, that "Henceforth 
space by itself and time by itself, are doomed to fade away into mere shadows, 
and only a kind of union of the two will preserve an independent reality." 
The importance of this is supported by the simplicity obtained by formulating 
the physical laws in four-dimensional space-time. In special relativity this 
simplicity is obtained by using the Minkowski space-time: 

ds^ = Tjpqdx'^dx'^ (1-3) 

where r]pq = diag{l, —1, {Throughout we use the convention that 
summation occurs over dummy indices, Latin indices take values from to 
3 and Greek indices values from 1 to 3, and take G — 1 and c — 1 units. 
The comma and semi-colon indicate ordinary and, covariant differentiation, 
respectively.) The presence of gravitation necessitates a generalization of the 
Minkowski space-time into the four dimensional Riemannian space-time: 

ds'^ = Qpqdx^dx^ (1.4) 

where Qah = Qahix^) is the metric tensor, symmetric in its indices, which 
characterizes the space-time completely. In the rest of this section we present 
important formulas/properties in Riemannian spaces required for our later 
work. 

A contravariant metric tensor g°-^ is defined as 

9""' = — , (1-5) 

where A"** is a cofactor of Qab, while g — d.ei{gpq). Prom this definition it is 
obvious that g"-^ is also symmetric in its indices and that: 

9bp9''' = 8t. (1.6) 

The relationship: 

dg = gg'"^ dg^^ = -ggpq dg^"^ (1.7) 

follows from (1.5) and (1.6). The metric tensor gab transforms under coordi- 
nates transformation — > as: 

dxP dx'^ 
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whereas its determinant g transforms under coordinates transformation as a 
scalar density of weight +2: 



dx 
dx 



9, 



[1.9) 



where 



dx 



dx 



is the Jacobian of the coordinates transformation. In general, a 



quantity which transforms like: 



dx 



dx' 



w 



dx"'' dx'^ 



dxP"'dx^"' « 



;i.io) 



is called a tensor density of weight W . Using dx"-' = ^^dx^ and (1.9) we can 

deduce that (—5') "'^^ 3?^;' " is an ordinary tensor if 3fJ^,'- - is a tensor density 
of weight W. Thus for a four-dimensional volume element d'^x the quantity 
y/—g d'^x is an invariant. An important tensor density of weight +1, called 
the Levi-Civita contravariant tensor density, is defined as: 



..abed 



+1, if abed is an even permutation of 0123 
—1, if abed is an odd permutation of 0123 
0, otherwise. 



(1.11) 



^abed jg totally auti-symmctric and has the useful property that its compo- 
nents are the same in all coordinate systems. Any totally skew-symmetric 
tensor of order 4 is proportional to this tensorial quantity. The covariant 
component Levi-Civita tensor density Eahcd may then be defined in terms of 
the contravariant component by the usual lowering of indices as: 



^abed — {—9)9ap9bq9er9ds£' 



pqrs 



(1.12) 



The tensor density Sabcd will be of weight —1. The above tensor densities 
can also be used to define the following totally antisymmetric unit tensors of 
rank four. The contravariant tensor e"'"^'^ is defined as: 



^abed 



-abed 



-9 



;i.i3) 



while the corresponding covariant tensor is defined by: 



^abed 



-9 ^abcd- 



:i.i4) 
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If Aij is an antisymmetric tensor, then using e"'^'^ and e'^^"^ we may define: 

*A'J = ^e^'^'Apg, (1.15) 

and 

*A'^ = ^e'^^'^A,,, (1.16) 

where *A^^ and *A'^^ are said to be a (i'ua/ pseudotensor and a dita/ tensor 
respectively to A^j. Using A*-' both and *Aij may be defined in a similar 
way in terms of Sabcd and Cabcd respectively. 

The motion of an infinitesimally small particle moving in a gravitational 
field is described by the geodesic equation: 

d'^x" „„ dxPdx'^ 



+ r" = fi 17) 

ds^ ^ P'' ds ds ' ^ ' 

which can be deduced from (1.4) by taking the deviation of the action inte- 
gral: 

7 = J ds, (1.18) 
where T^^ are the Christoffel symbols of the second kind given by: 

r?c = Ig"" i9bp,c + 9cp,b - 9bc,p) . (1.19) 
By contracting the pair of indices a and c in the above (1-19) we obtain 

= \9''9p,,b. (1.20) 

and combining this result with that of Eqn (1.7) we deduce that: 

1 dq 9 In -v/— Q 

2gdxb dx^ ■ ^ ' 

The identity that the g"-^.^ = gives us the following relationship: 

g-\^^-g<^VVl-gV^Vl,. (1.22) 

^Prom the definition of covariant differentiation we may deduce the following 
useful relation: 

Aa-b - A;a = ^a,6 " A,a- (1-23) 
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The Riemann curvature tensor is defined using the above Christoffel sym- 
bols as: 

pa -pa -pa , -pa -pp -pa -pp /i nA\ 

^ bed — ^ bd,c ^ bc,d "r pc^ bd pd^ bci \^-^^) 

which vanishes if and only if the space-time under consideration is fiat. In 
geodesic coordinate system (i.e. coordinate system where components of 
Christoffel symbols vanish at the pole of this special coordinate system) the 
Riemann curvature tensor takes the simpler form: 

-^Vd — ^bd,c ~ ^bc,d- (1-25) 

By lowering the upper index of the Riemann tensor (1.24) to Rabcd — 9apR^bcd 
then its symmetry properties become easy to study. Rabcd is anti-symmetric 
in each of the pair of indices ab and cd, and is symmetric under the exchange 
of these pairs of indices with each other. The cyclic sum of components of 
Rabcd formed by a permutation of any three indices is equal to zero. Several 
important tensors may be constructed using this tensor. By contraction of 
the Riemann curvature tensor we obtain the Ricci tensor defined by 

Rab = R^apb = 9'"^Rpaqb (1.26) 

— rp _ pp _L F9 FP _ ri rP n 07) 

~ ab,p ^ ap,b ^ ^ ab'- pq ap'- bq-> K'-'^' ) 

which is symmetric in its indices. The Ricci scalar curvature is: 

R^R% = gPm,,. (1.28) 
The Einstein tensor is defined in terms of the Ricci tensor as: 

Gab = Rab — -dab R- (1-29) 

Similarly the above Einstein's tensor is symmetric in its indices. The traceless 
Weyl conformal tensor Cabcd is again defined in terms of Riemann tensor as: 

Cabcd — Rabcd — ^ {OacRbd — QadRbc — QbcRad + QbdRac) 

+ ^ {9ac9bd - 9ad9bc) ■ (1-30) 
D 

The symmetry properties of Cabcd are similar to those of the Riemann cur- 
vature tensor. The importance of the Weyl tensor is in its invariance under 
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conformal mapping of two Riemannian spaces. The Bianchi identities are 
given by 

VaR',,, + V,i?V + Vci?V = 0- (1-31) 
These identities take a much simpler form if expressed in terms of the dual 
*Rpqrs to the Riemann tensor Rahcd 

VP*Rabcp^O. (1.32) 

By using the fact that the covariant derivatives of the metric tensor vanishes 
and contracting the Bianchi identities 1.31 we get the following important 
relationship: 

VpG/ = 0, (1.33) 
or the following equivalent identity: 

1 n D 

(1.34) 



1.2 Special relativity 

The idea of relativity can be traced back to Galileo who was the first to 
state clearly the concept of relative motion. His example of a boat which 
is in uniform motion illustrates that no experiment performed in a sealed 
cabin can be able to detect motion. According to the Galilean principle of 
relativity all laws of mechanics are the same in all inertial reference frames. 
The space a;" — > x"' and time t t' coordinates transformations between 
any two inertial frames, K and K\ may be written as: 

t' = t + to, (1.35) 

where constants M*^, v", and respectively represent rotations, uni- 
form motion, space and time translations of origins of reference systems. 
The so-called Galilean transformations (1.35) form a ten-parameters group 
of, namely three rotation Euler angles in the orthogonal matrix M, three 
components in each of and rf", and lastly the constant to- Newton's fun- 
damental equations of classical mechanics, based on the concepts of absolute 
space and absolute time, are invariant in all inertial reference frames under 
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Galilean coordinate transformations. If we consider three inertial coordi- 
nates frames fCi, and K-^^ then using the group property of Galilean 
transformations we obtain the following velocity addition law: 

<3=<2 + ^23, (1-36) 

where w^j, is velocity of Ki, relative to Ka- This velocity addition law (1.36) 
is vahd in Newtonian mechanics. However, Maxwell's theory of classical 
electrodynamics implied that light travels in vacuum at a constant speed 
c. Obviously, according to Galilean velocity addition law light could not 
have the same speed c with respect to arbitrary inertial frames. One can 
also verify that the wave equation for electrodynamics in free space is not 
invariant under Galilean transformations. The fact that Newton's equations 
of mechanics are invariant under Galilean transformations while the Maxwell 
equations are not, leads one to enquire whether it is possible to find some 
principle of relativity which holds for both mechanics and electrodynamics, 
but where 

• Newton's laws of mechanics are not correct, or 

• Maxwell's laws of electrodynamics are not correct? 

Experimental evidence indicating deviations from either the Newtonian or 
Maxwellian theory is required in order to decide whether the laws of me- 
chanics or electrodynamics need to be reformulated. 

The famous Michelson-Morley experiment showed that the velocity of 
light is the same for light travelling along the direction of the earth's orbital 
motion and transverse to it. Numerous attempts were made to explain the 
Michclson-Morlcy null result of the earth's motion through ether. Fitzgerald 
and Lorentz independently suggested that material bodies contract in the 
direction of the their motion by a factor ^1 — v"^ /c?. Lorentz explained the 
contraction hypothesis in terms of an electromagnetic model of matter, and 
introduced length transformation and the concept of a 'local time'. Poincare 
called for the development of a new mechanics to replace the Newtonian 
mechanics (see in Schroder [57]). He showed that the Maxwell equations in 
vacuo are invariant under the Lorentz transformations. He further showed 
that Lorentz transformations form a group [5 7]. The inhomogeneous Lorentz 
space-time x"- — > x°-' coordinate transformations , also called the Poincare 
group, are given by: 

x"' = A^x^' + d", (L37) 
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where constants d"- denote space and time translations, while represent 
rotations and uniform motion of the origins of inertial frames. The matrix A 
must satisfy the following condition: 



KKVpq = Vab, 



(1.38) 



with 



Vpq = V 



diag{l, —1, —1, — 1). 



(1.39) 



The so-called proper Lorentz group is obtained by imposing the following 
additional conditions: 



If there are no translations, i.e. d"' — 0, we get the homogeneous Lorentz 
group. 

The incompatibility of mechanics with electrodynamics was finally re- 
solved by Albert Einstein who proposed a replacement of the Galilean trans- 
formations by the Lorentz transformations. Einstein's theory of special rel- 
ativity (STR) is based on the following postulates: 

• The principle of relativity. The laws of physics are the same in all 
inertial reference frames. 

• The constancy of the speed of light. The speed of light in vacuum 
is the same for all inertial observers irrespective of the motion of the 



Following Einstein's approach, then the above Lorentz transformations (1.37) 
become a consequence of the constancy of the speed of light. As already men- 
tioned, in special relativity physical laws are formulated in four-dimensional 
Minkowski space-time. The interval Iab between two events A and B, with 
coordinates {Ia, xa, Ua, za) and {ts, xb, Ub, zb), given by: 



Iab = (4 - - (4 - 4)' - (4 - 4)' - (4 - 4)', (1-41) 



is an invariant under Lorentz transformations. (Note that = ct and that 
c = 1). Iab is said to be timelike if Iab > 0, spacelike if Iab < 0, and lightlike 
or null if Iab = 0. For a timehke interval of two events, it is possible to do 
a Lorentz transformation to an inertial frame where the two events occur at 
the same point, but there is no inertial frame where they occur at the same 
time. Therefore, for any two timelike events if > x^ then A is said to 



A° > 1, det{A) = 1. 



(1.40) 



source. 
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be in the absolute future of B, while if < then A is said to be in the 
absolute past of -B. In a spacelike interval for any two events it is always 
possible to do a Lorentz transformation to an inertial frame where the two 
events occur at the same time, but it is impossible to choose a frame where 
they occur at the same point. 

In Minkowski space-time the motion of a mass point is described by the 
world line = a;* (A), where A is parameter determining motion. Now, by 
observing that the proper time dr, given by: 

dr^ ^rjpqdx^dx", (1.42) 

is an invariant under Lorentz coordinate transformations,the velocity four- 
vector u'' and the acceleration four-vector are defined with respect to the 
proper time as: 

= f • 

and 

■■ du" /-, ..X 

= — . (1.44) 
or 

Newton's equations, in their original form, are no longer invariant under 
the Lorentz transformations, but using the above definitions the four-vector 
relativistic force of Newton's second law takes the form: 

^2 a 

r = rn^, (1.45) 
or in terms of the relativistic energy-momentum: 

p» = m^, (1.46) 
or 

the relativistic force may be expressed by: 

r^f. (1.47) 

where p"" is also a four-vector. /" is commonly referred to as the Minkowski 
force. The spatial components of the four- vector Minkowski force = (/°, f) 
are related to the Newtonian force F by the following: 

f = 7F, (1.48) 
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and its time component is given by 

f = jF- V, (1.49) 

where v and 7 := (1 — v^)~^^^. The time component of is the energy: 

= E = m7, (1.50) 

and the space components p° form the momentum vector: 

p = m7v. (1-51) 

Thus the four- vector p"' is called the energy-momentum vector. This energy- 
momentum vector is conserved in all inertial reference frames related by 
Lorentz transformations. For a particle with rest mass m > we have that: 

PaP" = m^ (1.52) 

which gives the following equation connecting the energy and momentum^: 

- p2 = m\ (1.54) 

Taking the ratio of equations (1.50) and (1.51) we get the following useful 
expression: 

I = V. (1.55) 

The Maxwell equations of electrodynamics in vacuum: 

V-E = 47rp, (1.56) 

VxB-— = 47rj, (1.57) 

V-B = 0, (1.58) 

VxE+— = 0, (1.59) 



^If wc now make use of the qiiantum mechanical correspondence between numbers 
and operators: pi— »— iV, E 1-^ —i-^ in (1.54) then one obtains the Klein- Gordon wave 
equation 

{D + m^)(t>{x) = 0, (1.53) 

for a scalar field 6 with mass m. 
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where E and B are respectively the electric and magnetic fields, are in- 
variant under Lorentz coordinate transformations. This may be illustrated 
by expressing the above equations in a covariant form in terms of an anti- 
symmetric electromagnetic field tensor Fab which may be defined as: 



?cib ._ 






Ex 


Ey 


E, 


—Ex 







-B 


-Ey 







Bx 


-E, 


By 


—Bx 






(1.60) 



Then using Eq. (1.60) the Maxwell equations (1.56, 1.57, 1.58, 1.59) may now 
be reduced to the following two tensorial equations: 



OF' 



ah 



= 0, 



(1.61) 
(1.62) 



where *F"'^ is the dual tensor of F""^, is the four-vector of current density 
with the components: 

r = (APv), (1.63) 

with p being the charge density, and v the velocity of charge. The equation 
of continuity for the current density: 



dx°- 



= 0, 



(1.64) 



now follows directly from Eq.(1.61). Eq.(1.62) can also be expressed as: 



Odd 

Fbc + TrrEca + TT-zFab = 0. 



dx" 



dx^ 



dx'^ 



:i.65) 



The relativistic electromagnetic force on a charged particle may be written 
in terms of the Maxwell tensor as: 



fa _ pa ^ 

where e is the charge parameter. 



;i.66) 
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1.3 The Energy- momentum tensor 

The energy-momentum of the electromagnetic field in free-space may be ob- 
tained by considering the action integral formed with the Lagrangian density: 

= (1-67) 

which gives the following contravariant components of the energy-momentum 
tensor: 

1 dA'P 1 

where A" is the four- vector of the potentials. This tensor is then symmetrized 
by adding the quantity and using the fact that in the absence of 

charges then Eq.(1.61) takes the form: 

which finally leads to the following symmetric energy-momentum tensor ex- 
pression: 

T"^ = J_(_i?»PF^ + -y^FpqPP'i). (1.69) 
47r 4 

This tensor (1.69) is gauge invariant and has a vanishing trace, T\ = and 
most importantly it satisfies: 

Tla = 0. (1.70) 
The fact that T\ is divergeless indicates that the four-vector: 



J a 



is conserved, where a is an arbitrary spacelike hypersurface. Now since the 
integral (1.71) is independent of the hypersurface cr, choosing the hypersur- 
face — constant, i.e. a three-dimensional space, then: 



may be identified as the total energy-momentum four- vector of the field. The 
angular momentum tensor: 
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defined in terms of T"^, satisfies: 

= 0, 

which shows that the total field angular momentum 

jab ^ f M'^'^^d^x, (1.74) 



is also conserved. 

The energy-momentum tensor (1.69) is only defined for charge-free fields. 
In the presence of charged objects then we should not only consider con- 
tributions from the electromagnetic field but also from the charged objects 
themselves. Thus the energy- momentum tensor of field together with charges 
is: 

rpab rpab i rpab 

J- — -t F "T -'M) 

where subscripts F and M indicate field and matter contributions, respec- 
tively. Now we should have that: 

^ (T^^ + T^^) = 0, (1.75) 



where Tp^ is given by (1.69). Now differentiating equation (1.69) and sim- 
plifying, and making use of the inhomogeneous Maxwell equations (1.61) 
together with (1.62) we get 



dx" 



For a system of non-interacting particles, the energy momentum tensor is 
given by 

d<i 

= (1.77) 

where /i defined as: 

=^ mn5(r - rn), (1.78) 
n 

may be termed the 'mass density' since it indicates the continuous mass dis- 
tribution in space. Noting that mass should be conserved for non-interacting 
particles, we have: 

d , ds. / „ X 

(/.-) = 0, (1.79) 



9x« ^ dt 
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and making use of Eq. (1.66) we get 



g^T^ - F'^J,. (1.80) 

Hence equations (1-76) and (1.80) show that equation (1.75) is indeed sat- 
isfied. Prom this discussion it is clear that electromagnetic fields are de- 
termined by the motion of charges, while on the other hand the motion of 
charges is determined by the fields, i.e. the two systems are interdependent 
on each other. 



1.4 The Einstein field equations 

In this section we give a derivation of the Einstein field equations for gravita- 
tion from a variational principle, but before that we give both the equivalence 
principle and the principle of general covariance which guided Einstein in his 
search for equations of gravitation. The equivalence principle is based on 
the equality of the gravitational mass and the inertial mass. Several experi- 
ments have been performed, starting from Galileo, Newton, Bessel, Eotvos, 
and many others, to investigate whether or not there is a difference between 
the gravitational and inertial masses of objects. The Eotvos torsion bal- 
ance experiment showed with a high degree of accuracy the equality of the 
inertial and gravitational mass. The principle of equivalence may be formu- 
lated as[75]: At every spacetime point in an arbitrary gravitational field it is 
possible to choose a "locally inertial coordinate system" such that, within a 
sufficiently small region of the point in question, the laws of nature take the 
same form as in a Minkowski spacetime , i.e. the gravitational field may be 
eliminated locally by the use of a freely falling coordinate system. 

The principle of general covariance may be stated in one of the following 
forms [9], namely 

1. All coordinate systems are equally good for stating the laws of physics. 

2. The equations that describe the laws of physics should have tensorial 
forms and be expressed in four-dimensional Riemannian spacetime. 

3. The equations describing the laws of physics should have the same form 
in all coordinate systems. 
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One considers the Ricci scalar curvature R — g^'^Rpq as our Lagrangian to 
derive the Einstein field equations. The action integral for the gravitational 
field is then given by: 

^Rd^x, (1.81) 



/ 



where the integration is taken over all space and over the time component 
x° between two given values. To encompass non-gravitational fields in a 
physical system we include another Lagrangian Lp for all the other fields, so 
that our expression for the action integral is given by: 



;i.82) 



where k is Einstein's gravitational constant. We require that the variation 
of the above be equal to zero, 

SI = 0. (1.83) 
The variation of the first part of I is: 

^ J V^Rd^x = J ^gP'i5Rpqd^x + J Rpq5{^gP'i) d^x. (1.84) 

By simpUfying and applying Gauss theorem it can then be shown that: 

j 5Rvq = 0. (1.85) 

Therefore 



S J V^Rd^x = I ^f^Rpq5gP''d^x + j R5^/^'"'d^x 

^ JV^ (Rp, - Idpg R ) SgP'i d'x. (1.86) 



Now the second part of action integral: 



-g Lpd X 



5g- 



pq 



dgPl dg'^S 
^T,,6gP'^d'x, 



d X 



;i.87) 
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where Tpq is the energy-momentum tensor defined as: 



T 



-g 



(1.88) 



Qgpq Qj.a 1^ QgP^ 

Now using results of (1.86) and (1.87) then the variation of (1.82) is: 

SI = J ^ (^Rpq - ^gpg R-KTpg^ SgP'' d^x, (1.89) 

and since 6g^'^ is an arbitrary variation this gives us the Einstein field equa- 
tions: 



(1.90) 



1.5 Maxwell equations in presence of gravi- 
tation 

In the absence of gravitation, the Maxwell field equations are expressed by 
(1.61) and (1.62). However, in the presence of gravitation these equations 
are generahzed by the following equations by replacing partial derivatives to 
covariant derivatives. Therefore one has now 



and 



The above two equations can easily be expressed as 



and 



d 



The continuity equation (1.64) now takes the form: 



(1.91) 
(1.92) 

(1.93) 
(1.94) 



= 



(1.95) 
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and the equation of motion of a charged particle in combined electromagnetic 
and gravitational fields, given as: 



m 



eF; 



ab 
b ' 



:i.96) 



is obtained from (1.66) by replacing du^-fdr with the intrinsic derivative 
Du"" /dr. Equations (1.91) and (1.92) arc the Maxwell equations in the pres- 
ence of gravitation. We now consider the Einstein field equations: 



pqi 



(1.97) 



in the presence of electromagnetic field, with the energy-momentum tensor 
Tpq for the electromagnetic field given by 



T 



pq 



-9 



-g Li 



d {d 



dg 



pq 



dx"- 



-9 Li 



The Lagrangian is given by: 



167r 



(1.98) 



(1.99) 



with the indices of Maxwell tensor components (see (1.92) for definition) 
Fab being raised or lowered by using metric components g"'^ or gab- Using 
the above, the energy-momentum tensor for electromagnetic field may be 
simplified to: 



ab 



1 

471 



FapF^ -|- ^gabFpqF^'^ 



(1.100) 



As this tensor is traceless, the Ricci scalar curvature R vanishes. Therefore 
the Einstein field equations in the presence of electromagnetic field are given 
by: 

Rab = f^Tab. (1.101) 



1.6 Klein- Gordon equation in curved space- 
time background 

In the absence of gravitation the Lagrangian density for a massive scalar field 
(j) is given by 

L=\ {v''M,q - • 
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The field equation obtained from above by using variation procedure is then 
given by 

dL d dL 

0, 



dxf d {dcp/dx^') 



which simphfies to 



□0 + m^(/) = 0, (1.102) 

where the differential operator □ := rj^'^d'^ / dx^dx'^ . This gives the familiar 
Klein-Gordon equation in a fiat space-time. To obtain the Klein-Gordon 
equation in curved space-time we invoke the principle of equivalence and 
general covariance. In the presence of gravitation the above equation (1.102) 
can be written as: 

p + mV = 0. (1.103) 
This is the Klein- Gordon equation in curved space-time. 



1.7 Conclusion 

In this introductory chapter we gave a list formulas which will be required in 
our later work. In the next chapter we will discuss the conservation laws and 
the dilemma associated with the definition of gravitational energy in general 
relativity. The lack of a generally accepted definition of gravitational energy 
has lead to doubts concerning energy localization in GR. The large num- 
ber of available pseudotensorial expressions used for computing energy and 
momentum distributions has even lead to suspicions that these nontensorial 
quantities would give different energy distributions for a given space-time. 
However, the pioneering work of Virbhadra on energy localization, on partic- 
ular space-time manifolds, has consistently shown this to be fallacious. In this 
work we are extending Virbhadra' s work by considering further space-times 
and showing that different energy-momentum complexes give the same en- 
ergy distribution in a given space-time. Hence energy-momentum complexes 
are useful expressions for computing energy distributions in GR. 
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Chapter 2 

Energy Localization 



2.1 Introduction 

The concept of total energy and momentum in asymptotically flat space-time 
is unanimously accepted; however, the localization of these physical quan- 
tities still remains an elusive problem when one includes the gravitational 
field. In the special theory of relativity the energy-momentum conservation 
laws of matter plus non-gravitational fields are given by: 

Tt,, = 0, (2.1) 

where T^^ denotes the symmetric energy-momentum tensor in an inertial 
frame, whereas general relativity leads to the following generalization of 
Eq.(2.1): 

= 0. (2.2) 

In this form Eq.(2.2) does not give rise to any^ integral conservation law 
whatsoever. In fact, if this equation is written as: 

with 



^However, Eq.(2.2) is the statement of the energy-momentum conservation laws in 
Special relativity in non-Cartesian and/or non-inertial coordinate systems. 
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then it is clear that the quantity Ki is not a general four-vector, therefore 
in a local system of inertia we can always make Ki to vanish in a given 
spacetime point and, in this case Eq.(2.3) simply reduces to Eq.(2.1). In 
general, ^ Q and for i = then Eq.(2.3) expresses the fact that matter 
energy is not conserved. 

Einstein formulated the conservation law in the form of a divergence 
to include contribution from gravitational field energy by introducing the 
energy-momentum pseudotensor t*^, so that: 

+ (2.5) 

The quantity t*^ is homogeneous quadratic in the first derivatives of the met- 
ric tensor and thus it is obviously not a tensor. With a suitable choice of a 
coordinates system can be made to vanish at a particular point. It can also 
be shown that if we form the integral J t^d^x in a flat spacetime using quasi- 
Cartesian coordinates,then its value is zero while if we transform to spherical 
coordinates the value of this integral is infinite (Bauer, 1918). Furthermore, 
it is possible to find a coordinate system for the Schwarzschild solution such 
that the pseudotensor vanishes everywhere outside the Schwarzschild radius 
(Schrodinger, 1918). Einstein ascribed these shortcomings to the coordinates 
used. However, the difficulties associated with Einstein's nontensorial quan- 
tities posed serious problems concerning the localizability of energy in general 
relativity. 

The problem of energy-momentum localization has been a subject of 
many research activities dating back to the very onset of the theory of 
general relativity but it still remains an open question. The numerous 
attempts aimed at finding a more suitable quantity for describing distri- 
bution of energy-momentum due to matter, non-gravitational and gravita- 
tional fields resulted in more energy- momentum complexes, notably those 
proposed by Landau and Lifshitz, Papapetrou, M0ller, and Weinberg. The 
physical meaning of these nontensorial (under general coordinate transforma- 
tions) complexes have been questioned by some researchers (see references in 
Chandrasekhar and Ferrari [15]). There are suspicions that different energy- 
momentum complexes could give different energy distributions in a given 
space-time. The problems associated with energy-momentum complexes re- 
sulted in some researchers even doubting the concept of energy-momentum 
localization in GR. According to Misner, Thorne and Wheeler [45] the energy 
is localizable only for spherical systems. However, Cooperstock and Sarra- 
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cino [22] refuted this viewpoint and stated that if the energy is locahzable in 
spherical systems then it is also locahzable for all systems. Bondi [6] wrote 
" In relativity a non-localizable form of energy is inadmissible, because any 
form of energy contributes to gravitation and so its location can in principle 
be found.'" It is rather unfortunate that the controversy surrounding energy 
localization which first appeared in electromagnetism,u;/iere there is an am- 
biguity concerning the choice of the Pointing vector, is also present in the 
most beautiful theory of general relativity (see in Feynmann, Lcighton and 
Sands[25]). The ambiguity in electromagnetism is not nearly as great a prob- 
lem as in general relativity because in the former, we are dealing with truly 
tensorial quantities. 

Over the past two decades considerable effort has been put in trying 
to define an alternative concept of energy, the so-called quasilocal energy. 
The idea in this case is to determine the effective energy of a source by 
measurements on a two-surface. These masses are obtained over a two-surface 
as opposed to an integral spanning over a three-surface of a local density as 
is the case for pseudocomplexes. A large number of definitions of quasi-local 
mass have been proposed, notable those by Hawkins, Penrose, and many 
others (see in Brown and York [8], Hayward[31]). Although, these quasi- 
local masses are conceptually very important (as Penrose emphasized) they 
still have serious problems. Bergqvist [5] furnished computations with seven 
different definitions of quasi-local masses for the Reissner-Nordstrom and 
Kerr space-times and came to the conclusion that no two of these definitions 
gave the same result. Moreover, the seminal quasi- local mass definition of 
Penrose is not adequate to handle even the Kerr metric (Beinsten and Tod 
[4]). On the contrary, several authors studied energy- momentum complexes 
and obtained stimulating results. The leading contributions of Virbhadra 
and his collaborators (Rosen, Parikh, Chamorro, and Aguirregabiria) have 
demonstrated with several examples that for a given spacetime, different 
energy-momentum complexes show a high degree of consistency in giving 
the same and acceptable energy and momentum distribution. In the rest of 
this chapter we present a brief introduction to each of the following : the 
Einstein, Landau and Lifshitz, M0ller, Papapetrou, and Weinberg energy- 
momentum complexes which are going to be used in our later work. 
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2.2 Einstein energy- moment una complex 

In order to arrive at Einstein's conservation laws for a system consisting 
of both matter and gravitational field we start with the gravitational field 
equations 



(2.6) 



and then using the contracted Bianchi identity {^W^ — |fl'*^-R).^ = 0, Eq. 
(2.2) becomes a consequence of the field equations. Eq. (2.4) can be written 
as 

(2.7) 



where g^^ = Using the field equations (2.6) we eliminate Tik from (2.7) 
and write Ki as 



16tt 



-99 



dL 



PI 



-.pi 



d 



dL 



-,pq 



dg 



where 



-9ti 



167r 



X^T _ ^pq 



and L is the Lagrangian density 



(2.8) 



(2.9) 



(2.10) 



Obviously i^*^ is a function of the metric tensor and its first derivatives. Now 
combining (2.8) with (2.3) we get the following equation expressing Einstein's 
conservation law: 

^"'^ 0, (2.11) 



dx^ 



where 



-9 {T.' + U') 



(2.12) 
(2.13) 
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is the total energy-momentum complex for the combined matter plus grav- 
itational field, while $5^'^ = ^—gT.^ ^ and v^^ = \/—gt.i^. By introducing a 
local system of inertia, the "gravitational" part v^^ can always be reduced 
to zero for any given space-time point. In general Q'j'^ is a function of matter 
and gravitational tensor, and hence the division of 9^^ into "matter" part 
and "gravitational" part is highly arbitrary. The matter part may even be 
eliminated entirely from (2.12) and the 9^^ expressed only as a function of 
the metric tensor together with its first and second derivatives, as: 

e=^, (2.14) 

with S^'^^ given in Tolman[61] as 

M0ller [46] suggested a more useful expression for 9^^. A quantity 9^^ which 
satisfies Eq. (2.11) identically, must be writable in the form: 

, ^ (2.16) 

167r dxP ' ^ ' 

where h^'^^ = —h/'^ and h-'^^ is a function of the metric tensor and its first 
derivatives. It is easy to verify that this is the case for 

h^^ = ^ [{-g) (^'V"^ - ^' V"^)] ,^ . (2.17) 

If the physical system under consideration is such that we can introduce 
quasi- Cartesian coordinates for which the gik converge sufficiently rapidly 
towards the constant values rjik where: 

riik^diag(l, -1,-1,-1) (2.18) 

then it follows from Eq. (2.11) that the quantities: 

Pi^ [ [ [ 9i^dx^dx^dx^ (2.19) 
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are constant in time, provided that 9^^ are everywhere regular. The integral 
in Eq. (2.19) is extended over all space for x° = const. Further, Gauss's 
theorem furnishes 

Pi-^J J K'^'i^adS (2.20) 

where /x^ = 7- is the outward unit normal vector over an infinitesimal surface 

element dS 

The main problem with interpreting the integrand 9^ in Pq as the energy 
density is that it does not behave like a three-scalar density under purely spa- 
tial transformations. It can be shown that if we form the integral J t^d^x 
then its value is zero in a flat spacetime using quasi-Cartesian coordinates, 
while if we transform to spherical coordinates the value of this integral is 
infinite [46]. Furthermore, Schrodinger [46] showed that there exists a co- 
ordinate system for the Schwarzschild solution such that the pseudotensor 
vanishes everywhere outside the Schwarzschild radius. Many other prominent 
scientists, including Weyl, Pauli, and Eddington, questioned the nontenso- 
rial nature of t^.* because with a suitable choice of a coordinate system it 
can be made to vanish at any point in spacetime, (for details see in Chan- 
drasekhar and Ferrari [15]). Einstein (see in Goldberg [28]) pointed out that 
these effects were artifacts of the coordinates used and that they were not 
related to the physical system used. Einstein showed that for a spacetime 
that approaches the Minkowski spacetime at spatial infinity then the energy- 
momentum Pi transforms as a four-vector under all linear transformations. 
Consider the following example of a closed system at rest whose coordinates 
x"- are chosen so that at large distances the line element is given by 

ds^ = (1 - a/r) {dx'^y - (1 + a/r) [[dx^f + {dx^f + (dx^)') (2.21) 

with = {{x'')^ + {x^)^ + (,T''^)^j, and where a = 2Mq could be regarded as 

a constant connected to the total Newtonian gravitational mass Mq. Using 
Eq. (2.20) the total energy-momentum components in this system will be 
given by 

P, = = -51M,. (2.22) 

Now using a Lorentz transformation of the form 

^ x' + vx^ ^ x^ + vx' ^2' ^^2 ^3'^^3_ ^2.23) 
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the total energy-momentum components Pj/ in the new coordinates system 
x""' have values 



Pv = { , , , , 0, ^ 2.24 

I ^r^' ^r^' 'J ^ ^ 

in other words, in an inertial system the total energy-momentum components 
have the same values as the components of the four-momentum of a particle 
of proper mass Mq moving with velocity v along the x-axis. Another im- 
portant result given by Einstein is that any two systems of quasi-Cartesian 
coordinates S and Si which coincide at spatial infinity, but differ arbitrar- 
ily elsewhere will have P* = P*'. Although Einstein was able to show that 
the energy-momentum pseudo-complex 9^^ provides satisfactory expressions 
for the total energy and momentum of closed system in the form of three- 
dimensional integrals (2.19), to get meaningful values for these integrals one 
is restricted to the use of quasi-Cartesian coordinates. 

An alternative form of Einstein's pseudo-complex, which we found useful 
in some of our calculations, is given by Tolman [61] as: 



1 

,m 



(2.25) 



where L is the Lagrangian given by (2.10) whereas q"'^ :— ^/^g while 

b 



'^^g°'^ c 1 that: 



^--Kb + \^lK^\^t^lr (2-26) 

Another useful expression, from Tolman [61], for obtaining energy for a static 
or quasi-static system using quasi-Cartesian type of coordinates is: 

E = j j j - - c^^' - ^3^) dx^dxV. (2.27) 

The main advantage derived from using this expression is that it can be 
evaluated by integrating only over the region actually occupied by matter or 
electromagnetic energy, since Q'^'' vanishes in empty space. 
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2.3 Landau-Lifshitz energy- momentum com- 
plex 

One of the main objections to Einstein's energy- momentum complex was 
that it is not even symmetric in its indices, so cannot be used to define con- 
servation laws of angular momentum. In this section we discuss an energy- 
momentum complex which satisfies this requirement. In deriving the con- 
served total four-momentum for a gravitational field plus matter and all 
non-gravitational fields. Landau and Lifshitz [40] introduced a geodesic co- 
ordinate system at some particular point in spacetime in which all the first 
derivatives of the metric tensor gn- vanish. Then at this point Eq. (2.2) can 
be reduced into a form: 

1^ = 0. (2.28) 

similar to Eq. (2.1). Further, it can be shown that quantities T^^, which 
satisfy Eq. (2.28) identically, can be expressed in terms of the following 

= (2.29) 

where the quantities S'^^^ are antisymmetric in their last two indices k and I. 
At the point under consideration the Ricci tensor may be written as 

{9pr,qs + 9qs,pr 9pq,rs 9rs,pq ) , (2.30) 

since the Christoffel symbols vanish. Now using the gravitational field equa- 
tions (2.6) we may deduce that the energy-momentum tensor T''^ can indeed 
be expressed as 

where the expression inside curly brackets can be associated with S'^^K Defin- 
ing the quantities: 

h'"' = [(-g) {g^'g'"^ - g'V^)] , (2.32) 

then obviously /i**^' = — /i'"^and since all the first derivatives of the metric 
tensor vanish, Eq. (2.29) may be written as: 

fihikl 

-Q^-i-g)T'' = 0. (2.33) 
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Eq. (2.33) will only hold in spacetime at some particular point of a special 

coordinates system in which all the first derivatives of the metric tensor go. 
vanish. In an arbitrary coordinate system, the difference dh^^^ /dx^ — (— (?) T**^ 
will no longer be zero. We denote this difference by: {—g) f'^, and thus in 
general Eq. (2.33) will be of the form: 

(-^)(T^'= + f'=) = ^, (2.34) 

where the quantities t are symmetric in their indices since both T**^ and 
Qjfki jQx^ are symmetric in the indices i and k. It is obvious that is not 
a tensor quantity. T^^ can be eliminated from Eq. (2.34) by making use of 
the Einstein gravitational field equations , the above equation may then be 
written as 

Ak ^ ( jyik "^^i^D 1 I uikl 



(-9) = ^ (R'' - ^g'^Rj + h''' ,1- (2.35) 

Now using the expression of /i*'^' in Eq. (2.32) and that of the Ricci tensor, 
and after a lot of simplifications the expression of reduces to the following 

' y y V pa qr ' qr pa ra pq pq ra) 

I rJ^Pri'p^ ( r*^ -I- r*^ r*^ r*^ ^ 

' y y V pa qr ' qr pa ra pq pq ra) 

I kp qr / pi pa i pi pa pi pa pi pa \ 

' y y \^ pa qr ' qr pa ra pq pq ra) 

+ gp^g''" (r -r r*^ )| . 

I y y V pr qs pq rsj } 

Also, since /i*'^' is antisymmetric in indices k and /, it follows from equation 
(2.34) that 

-^[i-9){T'' + f')]=0, (2.36) 
which means that there is a conservation law for the quantities 

P' = J i-y) {T'^ + t'^) dSk- (2.37) 

where the integration may be taken over any infinite hypersurface includ- 
ing all of three dimensional space. In the absence of gravitation, in quasi- 
Cartesian coordinates system, the set of quantities t**^ vanishes and P* re- 
duces to J {—g) T'^^dSk which is the four-momentum of the physical system 
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without gravitation. Therefore P* in (2.37) is identified with the total four- 
momentum of the whole physical system including gravitation. So we refer 
to as the energy-momentum pseudo-tensor and to 

V'' = {-g) {T^ + f^) (2.38) 

as the energy-momentum complex. Now choosing the hypersurface — 
const, then can be written in the form of a three dimensional space 
integral 

P^= j j j L^dx^dx^dx^. (2.39) 

Hence we might interpret the quantity L^^ as representing the energy density 
of the whole physical system including gravitation, and interpret the quantity 
L^*^ as representing the components of the total momentum density. 

Unlike the Einstein energy-momentum complex 6f, the main advantage 
with the Landau-Lifshitz energy-momentum complex L*'^ is that it is sym- 
metric with respect to its indices, and therefore it can be used to define a 
conservation law for the angular momentum. We define it as 




= J {x'V^^ -x^L"'')dSm. (2.40) 

By using an argument similar to one used by Einstein, it can be shown that: 

1. For asymptotically fiat spacetime the quantities P' are constant in time. 

2. For any two systems of quasi-Cartesian coordinates S and St which 
coincide at spatial infinity, but differ arbitrarily elsewhere we have P' = 

3. P* transforms like contravariant components of a four- vector under all 
linear transformations, including Lorentz transformations. Therefore 
using (2.34) the energy-momentum densities (2.37) may be written as 

P'- j h'^'^^r^dSk = \j ih^'"',mdSk - h''"%dSm) (2.41) 
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so that the above integral can then be written as an integral over an 
ordinary surface giving 



where d*nkm is the normal to the surface element related to tangen- 
tial element dn^^^ by d*nik = ^eikmgdn^'^. Choosing the hypersurface 
= const for the surface of integration in (2.37) then the surface of 
integration in (2.42) becomes an ordinary space, thus we obtain 



where dum — d*nom is a three-dimensional element of an ordinary space. 
Similarly, an analogous formula for angular momentum is given by 



2.4 M0ller energy- momentum complex 

M0ller [46] argued that although the Einstein energy-momentum complex 
provides useful expressions for the total energy and momentum of closed 
physical systems, the singling out of quasi-Cartesian coordinates is somehow 
unsatisfactory from the general relativity viewpoint. Most of the criticism of 
Einstein's prescription centred around the nontensorial nature of the quan- 
tity t^^. A mere change of a coordinates system from quasi-Cartesian into 
spherical polar coordinates creates energy in vacuum. So M0ller searched 
for an expression of energy and momentum which is not dependent on any 
particular coordinates system. If 6-'^ is Einstein's energy-momentum com- 
plex and Sj^'' is another quantity with an identically vanishing divergence, 
then their sum 6- + S^'^ will also vanish identically. Therefore the energy- 
momentum complex is not uniquely determined by the condition that its 
divergence vanishes. M0llcr [46] exploited this freedom by searching for a 
quantity S^^ that can be added to 9^'' so that it transforms as tensor for spa- 
tial transformations. In order to retain the satisfactory features of Einstein's 
theory S^ ^ had to be chosen in such a way that it was form invariant function 
which depends on the metric tensor and on its first and second derivatives. 
Under linear transformations it had to behave like a tensor density satisfying 
the following conditions: 




(2.42) 




(2.43) 




,i I, fcOm 
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1. — identically, therefore it must be expressible in terms of ^j'^^p 
where ^'^'^^ = — ^j^'^is an affine tensor of rank 3. 

2. J Sj^^d^x — over total three-space for a closed system if we use quasi- 
Cartesian coordinates. 

3. 9q'' + Sq'' behaves like a four- vector density under all transformations 
of the type: 

x'^'^r'ix^). (2.44) 

Thus for 

Q,'' = 9,' + S,', (2.45) 

condition (2) implies that 

%°dx^dx^dx^ ^ J J J Oi'^dx^dx'^dx^, (2.46) 

for a closed physical system. In order to find an S^'' satisfying the above 
conditions M0ller [46] first investigated transformation properties of un- 
der arbitrary infinitesimal transformations of the type Eq. (2.44) so as to 
establish the deviation of the variation of 6q^ from a scalar density. Following 
this procedure he finally arrived at: 

where 

= 2h^' - S^hf' + 5lhf (2.48) 

dgip dgig\ . 




- v^^^-^J/V^ (2.49) 

Because of the antisymmetry, Xi'^ — ~Xi^-, it now follows that: 

ac^ k 

"^^^ 0. (2.50) 



^k 



X'^ 

M0ller was able to show that the quantities: 




%'^dx^dx'dx% (2.51) 
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for a closed system at rest in quasi-Cartesian coordinates coincide with the 
corresponding ones for the Einstein pseudo-complex, and hence equation 
(2.46) is satisfied and then following an argument similar to Einstein's he then 
deduced that Pi transform like four-vectors under Lorentz transformations, 
and thus the integrals have the same values as the corresponding integrals of 
Einstein in all cases where the latter are meaningful at all. Finally, he showed 
that $5q transforms like a four-vector under the transformations of the form 
(2.44). Using Gauss's theorem the total energy- momentum components are 
given by 

= /x.°>ac?^, (2.52) 

where /x^ is the outward unit normal vector over an infinitesimal surface 
element dS 

M0ller's coordinate independent prescription appeared to have finally 
solved the problem of energy localization until three years later when M0ller 
[47] performed a Lorentz transformation for a closed system at rest using 
the line element (2.21). As shown above (2.24), the Einstein pseudo-complex 
exhibited the correct transformation properties. M0ller [47] first obtained 
the total energy-momentum components Pj such that 

Pi = -5°Mo, (2.53) 

in the coordinates system x", using Eq. (2.51), which agrees with those 
obtained in Eq. (2.22) using Einstein's prescription. Now using a Lorentz 
transformation, Eq. (2.23), to obtain the total energy-momentum compo- 
nents Pi' in the new coordinates system he got the following values 

Mo A , 2 2\ 5 Mqv 



which did not seem to give expected transformation properties under the 
Lorentz transformation, and hence showing that Pj does not transform like 
a four- vector under Lorentz transformations. After a critical analysis of 
M0ller's result, Kovacs [39] claimed to have found a defect in M0ller's cal- 
culation. However, Novotny [50] has shown that M0ller[47] was right in 
concluding that Pj does not transform like a four-vector under Lorentz trans- 
formations. Lessner [41] finally showed that the problem lies with the inter- 
pretation of the result from a special relativistic point of view instead of 
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a general relativistic point of view. According to Lessner [41]: The energy- 
momentum four-vector can only transform according to special relativity only 
if it is transformed to a reference system with an everywhere constant velocity. 
This cannot be achieved by a global Lorentz transformation. He concludes by 
stating that M0ller's energy-momentum complex is a powerful expression of 
energy and momentum in general relativity. 



2.5 Papapetrou energy- momentum complex 

Amongst the five energy-momentum complexes under discussion, the Pa- 
papetrou energy-momentum complex is the least known and as a result it 
has been rediscovered several times, first in 1948 by Papapetrou. Gupta, 
a high energy physicist well-known for the Gupta-Bluer formalism in quan- 
tizing the electromagnetic field, in 1954 also obtained this complex using a 
slightly different method (sec Gupta [29]). This is the reason Misner refers 
to this complex as the Papapctrou-Gupta energy-momentum complex. In 
1994, a renowned particle physicist, R. Jackiw of MIT - USA, with his col- 
laborators re-obtained the same energy-momentum complex using the same 
method used by Papapetrou (see Bak, Cangemi, and Jackiw [2].) This fact 
was pointed out to them by Virbhadra and then they sent an errata to PRD. 
Papapetrou [51] formulated this conservation law of general relativity by ex- 
plicitly introducing in calculations and in the final formulae the fiat-space 
metric tensor r)ik- The formula was obtained following the generalized Be- 
hnfante method. First, from Eq. (2.18) Einstein's pseudocomplex may be 
written as: 



SndxP 

where 

^ fern 



-Q^ — — + -6- 



dg^^^ ' 2"- c}/VJ' ^^-^^^ 

and L is the Lagrangian given by (2.10), fp'^^and g^^^ are the same as defined 
previously. Now to symmetrize the above total energy-momentum complex 
using Belinfante method we start by assuming the existence of a quantity 
Jl*'^ so that Jl*'^ = Q*^' which differs from rf^Oj^ only by a divergence: 



V'^C + —B'^P (2.57) 
' p dxP ^ ' 
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so that its divergence vanishes: 

d 



Vl'^ = 0. (2.58) 



The above equation will only be satisfied if i?*^' is antisymmetric in its last 
two indices, i.e. if B'^^^ — —B"^^^. By making use of the Belinfante method 
^iki j^g^y |-jg expressed by the relation: 

^iki ^ _ 1 (^jiki _^ jUk _^ jiki^ ^2.59) 

where f^^^ is the spin density of the field given by: 

1 f)T 

f'"' - ^^y^^"" - ^"''^') ^ (2.60) 
which, using Eq. (2.56), may also be written as 

= ^ (^/^'^ - Ri'r^'') > (2-61) 
Now using this in Eq. (2.57) and (2.59), we get: 

" li^^ {RI' + Ra ") + {Rr + RD - r,"" {R:' + Ra"')] 

^ ^ (2.62) 
then using (2.26) in (2.56) we get 



b^cp 



Ra' + R:' = (s^aQ"' - Is'aQ'' " 

which is used to simplify (2.62) to 
where 

Note that the quantities jV*'^"^ are symmetric with respect to the first two 
indices i and k. The energy-momentum complex Q**^ of Papapetrou satis- 
fies the local conservation laws (2.58). This locally conserved quantity O**^ 
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contains contributions from the matter, non-gravitational and gravitational 
fields. and ^2"° are the energy and momentum (energy current) density 
components. The energy and momentum components are given by 

P'^ J J J n'^dx^dx^dx^ (2.65) 

and for the time-independent metrics Gauss's theorem furnishes 



where Ua is the outward unit normal vector over an infinitesimal surface 
element dS. 

The energy-momentum density Q'*^ is symmetric with respect to the two 
indices i and k, therefore it can be used to define angular momentum density 



2.6 Weinberg energy- momentum complex 

Weinberg [75] obtained an energy-momentum complex by considering a quasi- 
Minkowskian coordinate system^. In this quasi-Minkowskian coordinate sys- 
tem Qab may be viewed as the sum of the components of the Minkowski metric 
r]ab and that part hab which vanishes at infinity, i.e. 

gab = nab + Kb- (2.66) 

hob is only assumed to vanish at infinity but may take arbitrarily large values 
elsewhere. Using this notation the Ricci tensor may then be written as: 

Rab^R^\,-rR^\,-rO{h% (2.67) 
where the linear part R^^\i, in hab is given by: 

2 V dx^dx^ dxPdx'' dxPdx" dxPdxr ' ' 



^We refer to a coordinate system as being quasi-Minkowskian if the metric gab ap- 
proaches the Minkowski metric rjab far away from a given finite material system. 
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and the second-order part R^'^\^ in hah is given by: 



2 



pq 



aq 



+ 



ab 



dx^dx"- dx^dxP dx'^dx"- dx'^dx^ 



+ 



dh 



ab 



dx^ 
dx„ 



+ 



dx"- 
dhPi 



dxg 



dxP dxi 



dx^ dxr, 



dh 



pb 



+ 



dh 



pq 



dh 



qb 



dxi dx^ dxP 



In the above expressions of R'^^\i, and indices on hik and d/dx^ are 

raised and lowered with ry's, whereas indices on true tensors such as Rik are 
raised and lowered with ^f's as usual. Using the above notation, Einstein's 
field equations may now be written as: 



R^^\b - \riabR^^^\ = -StT [Tab + tab] , 



where 



1 

Stt 



Rah-\gabR% - R^'\,-\r)abR^'^% 



Note that the left hand side of Eq. (2.69) may also be written as 



where 



dh*^ dh"- 



dxi 



dXr. 



dh"' 
dx"' 



Qj^ap 

dx"' ' 



ik 



dh 



ik 



dXr. 



QJ^pk 

dxi 



(2.69) 
(2.70) 

(2.71) 

(2.72) 



is antisymmetric in its first two indices i and p, that is, uj^'^ — —uj'^^^ so that 



it now follows that: 



d 

dxP 



2 ' 



0. 



(2.73) 



After analyzing Eq.(2.69) and comparing its form with that of the wave 
equation of a field spin of 2, and further noting that, since the quantities 
R^^\ij obey the linearized Bianchi identities Eq.(2.73) Weinberg concluded 



that the quantity W given by: 



W 



ab 



V^'^V^lTpq + t 



pqi : 



(2.74) 
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which by Eqs. (2.69) and (2.73) satisfies the following local conservation law: 



d 

= 0, (2.75) 



dxP 



may be interpreted as consisting of the total energy-momentum pseudocom- 
plex of matter and gravitation, with tik indicating the energy-momentum 
pseudotensor of gravitation. Therefore for any finite system of volume V 
bounded by the surface S, we have: 

4- I W'Sx = - I W'riadS, (2.76) 
dt Jv Js 

where n is the unit outward normal to the surface. Analogously, the quanti- 
ties P* given by: 

W^'dx^dx^dx"^, (2.77) 




may therefore be interpreted as representing the total energy-momentum 
components of the system including matter, electromagnetism, and gravita- 
tion. 

Since hah — > as r — > oo, the energy-momentum tensor of matter plus 
non-gravitational field Tab also vanishes at infinity. Prom 



hab = O 



hab,c = ^ I ^ 



1 
r 
1 
r 

hab,cd = 0(^^J, (2.78) 

then using Eqs. (2.67,2.68,2.70) and noting that the quantity tab is of the 
second order in /i, it follows that 

ti'^ = ^ (^) ' (2.79) 

approaches zero at infinity. This shows that the source term on the left hand 
side of (2.69) is effectively confined to a finite region. Therefore the quantities 
in (2.77) that give the total energy and momentum converge. 
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Weinberg further justifies fiis choice of the quantities by showing that 

these quantities arc four-vectors and are additive. We illustrate below that 
arc invariant under any transformation that reduces to an identity trans- 
formation at infinity 

1. First consider the transformation of the form 

x"' = + r{x) (2.80) 

where /"(.x) ^ as r — oo. Then to first order in both and hab-, 
the metric tensor g°- ^ in the new coordinate system x'^ will then be 
given by 

ga'b' ^ ^ab _ ^a'b' (2.8I) 

where 

^a'b' = h^b_di^_df_ ^2.82) 

OXb OXa 

since as r — cxd both and hab are small. This change in coordinate 
transformation will lead to the following change in quantity uj^^'^ defined 
by (2.72) 

= P)«^^^„ (2.83) 

where 



2 [ dxi dxb dxa dxa dxb dxi 

(2.84) 

Let us express P'^ in terms of u;^^'^ as 

P*^ = +^ / / / "^^^^P^^^d^^^^^ (2-85) 
which, using Gauss's theorem, gives 

jfc , ^ f f , .aOk , 



P j ^ ' l^adS (2.86) 

where = ^ is the outward unit normal vector over an infinitesimal 
surface element dS = r'^ sin 6 dOd^p. Now noting that is totally 

antisymmetric with respect to its first three indices a, b and i, then the 



40 



above change Ao;^**^ will result in the following change in the surface 
integral 



AP" = +^ j j D^^^'^^ii^dS (2.87) 



87r 

+ 



87r 

which, using Gauss's theorem, gives 

AP^ = +— [ [ [ D^^^'^^f^fXadx^dx^dx'^ = 0, (2.88) 



thus showing that is invariant under any transformation that reduces 
to an identity transformation at infinity. An important consequence of 
this result is that transforms as a four- vector under any transfor- 
mation that leaves the Minkowski metric rjab at infinity unchanged, 
because any such transformation can be expressed as the product of a 
Lorentz transformation under which transforms as a four- vector. 

2. We now show the additive property of P*. Dividing the matter in our 
system into distant subsystems we can approximate the gravita- 
tional field hab as the sum of /i^^'* 's that would be produced by each 
subsystem acting alone. Thus, from the above calculation of P*, it fol- 
lows that the total energy and momentum are equal to the sum of the 
values P^^^-j for each subsystem alone. 

Now from the above, it has been shown that P* is conserved, is a four-vector, 
and is additive. In addition to these properties, the total energy-momentum 
complex W^'' is conserved and symmetric in its indices, therefore we can use 
it to define a conservation law for angular momentum 

M^i" = (2.89) 

where 

SO that M^'^ and M"**^ can be taken as representing the density and flux of 
a total angular momentum 

J^*^ = / / / M^'^dx^dx^dx^ 
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which is constant if M"**^ vanishes on the surface of the volume of integration. 

Both MP''' and J^^ are antisymmetric with respect to indices i and k. As 
above, the total angular momentum complex can be written in terms of a;^**^ 
as 

f f I {x'uj^''\-x''uj'''^\) dx'dx^dx'^ (2.90) 



which, using Gauss's theorem, gives 

J J { dx- dt 

Oh ^ 

+ + ^OpSa^ - hoaSp^ > l^adS. (2.91) 

We only give the physically meaningful components of J**^ which are the three 
independent space-space components = J^^ , = J^^ and = J^^. In 
order to calculate the total momentum, energy, and angular momentum of 
an arbitrary system, one only needs to know the asymptotic behavior of h^b 
at great distances. Though the quantities tab- W"'^ and M^*'^ are not tensors, 
they are at least Lorentz covariant. Thus for a closed system P* and J*^ will 
not only be constant but Lorentz-covariant. 



2.7 Present study of energy localization 

Rosen and Virbhadra [56] investigated the energy and momentum of the 
Einstein-Rosen metric using the Einstein energy-momentum complex. The 
Einstein- Rosen metric is a non-static vacuum solution of Einstein's field equa- 
tions that describes the gravitational field of cylindrical gravitational waves 
given in cylindrical polar coordinates (p, (p ,z) by the line element 

ds^ = e^^-^^'idt' - dp') - e-'^'pW - e'^'dz^ (2.92) 

where 7 = 7 (p, t) and ^ — ^ {p, t) and 

7,t = 2p^',p*,t, 

7,t = p{^l + ^l). 
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These authors carried out their calculations in quasi-Cartesian coordinates, 
and reported that the energy and momentum density components are non- 
vanishing and reasonable. Rosen (see in [54]) had earlier computed, in cylin- 
drical polar coordinates, the energy and momentum components in this met- 
ric using the energy-momentum complexes of Einstein and Landau and Lif- 
shitz. For both prescriptions the energy and momentum density components 
vanished. Initially, the vanishing of these components seemed to confirm 
Scheidegger's conjecture that a physical system cannot radiate gravitational 
energy. However, two years later Rosen [54] realized his mistake and recal- 
culated energy and momentum density components in quasi-Cartesian coor- 
dinates and found finite and reasonable results, which were later reported 
by himself and Virbhadra [56]. Virbhadra [70] showed that the energy- 
momentum complexes of Einstein and Landau and Lifshitz give the same 
energy and momentum densities when calculations are carried out in quasi- 
Cartesian coordinates. The energy density of the cylindrical gravitational 
waves was found to be finite and positive definite. The momentum density 
components was found to reflect the symmetry of the spacetime. 

In a recent paper Virbhadra [74] showed that the energy-momentum com- 
plexes of Einstein, Landau and Lifshitz, Papapetrou, and Weinberg, and the 
Penrose quasi-local definition give the same result for a general nonstatic 
spherically symmetric metric of the Kerr-Schild class. The well-known space- 
times of the Kerr-Schild class are for example the Schwarzschild, Reissner- 
Nordstrom, Kerr, Kcrr-Newman, Vaidya, Dybney et ai, Kinnersley, Bonnor- 
Vaidya and Vaidya-Patel. These spacetimes are defined in terms of the fol- 
lowing metrics: 

9ab = Vab - Hkah, (2.93) 

where rjab is the Minkowski metric, if is a scalar field, and ka is a null, 
geodesic and shear free vector field in the Minkowski spacetime. Each of 
these are given by 

r]''U:,kg = 0, (2.94) 
V^'kpK = 0. (2-95) 
iK,, + k„p)KrV'''-KY = 0. (2.96) 

The vector ka of the Kerr-Schild class metric Qab remains null, geodesic and 
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shear free with the metric Qab ■ Eqs. (2.94)-(2.96) lead to 



g^'^kpkq — 0, 

9^'^ki-pkq = 0, 

{kp-^q + kq-^p) k^.j.g'^ — (^"^p) — 0. 
Aguirregabiria et. al.[l] obtained the following results 

- I^^-^^'^^ (2-97) 

where 

^abcd _ ^ (^"fefc'^fc'i + r]^k''k^ - r]'^k^k'^ - rf'^k'^k") 

therefore the energy- momentum complexes of Einstein 9^, Landau and Lif- 
shitz U^, Papapetrou Q**^, and Weinberg W^^ 'coincide' for any Kerr-Schild 
class metric. Only the null conditions of equations (2.94) - (2.96) was used 
to obtain the above results in terms of the scalar function H and the vector 
ka for the Landau and Lifshitz, Papapetrou, and Weinberg complexes was 
used, while for the Einstein complex the null as well as geodesic conditions 
were used. These energy-momentum complexes 'coincide' for a class of solu- 
tions more general than the Kerr-Schild class since the shear-free conditions 
were not required to obtain the above equations. The energy and momentum 
components are given by 



IGtt 




,9 



Since the energy-momentum complexes of Landau and Lifshitz, Papapetrou, 
and Weinberg are symmetric in their indices the corresponding spatial com- 
ponents of angular momentum are defined as 

2.8 Conclusion 

In this chapter we elaborated on the problem of energy localization in gen- 
eral relativity, a concept which still remains a puzzle. There have been many 
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attempts at finding an appropriate method for obtaining energy-momentum 
distribution in a curved space-time, which resulted in various energy-momentum 
complexes. These complexes are restricted to the use particular coordinates. 
The problem associated with energy-momentum complexes resulted in some 
researchers doubting the concept of energy localization. However, the lead- 
ing contributions of Virbhadra and his collaborators have demonstrated with 
several examples that some of these complexes consistently give the same and 
acceptable energy and momentum distribution for a particular spacetime. We 
elaborated on only a few energy-momentum complexes used in our work. We 
also highlighted alternative attempts aimed at finding the alternative con- 
cept of quasilocal energy-momentum. The coordinate-independent quasilocal 
mass definitions are important conceptually; however, serious problems have 
been found with these. Chang, Nester and Chen[14] have also shown that the 
Einstein, Landau and Lifshitz, M0ller, Papapetrou, and Weinberg energy- 
momentum complexes may each be associated with a legitimate Hamilto- 
nian boundary term, and because quasilocal energy-momentum are obtain- 
able from a Hamiltonian then each of these complexes may also said to be 
quasilocal. 
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Chapter 3 

Energy Distribution of Charged 
Dilaton Black Hole Spacetime 

3.1 Charged dilaton black hole spacetime 



Properties of the well-known Reissner-Norstrom black holes have been stud- 
ied extensively. If one couples the dilaton field to the Maxwell field of the 
black holes many of these properties are known to change (sec in Holzhey 
and Wilczek [32]; Horne and Horowitz [33]). Charged dilaton black holes 
have been a subject of study in many recent investigations [32]- [34]. Garfin- 
kle, Horowitz and Strominger (GHS) [27] considered the action 



where R is the Ricci scalar, $ is the dilaton field, = F^bF"-'^ where Fab 
is the electromagnetic field tensor, and 7 is an arbitrary parameter which 
governs the strength of the coupling between the dilaton and the Maxwell 
fields. We shall consider only nonnegative values of 7 since a change in the 
sign of the parameter 7 will have the same effect as a change in the sign of 
the dilaton field $. The action (3.1) reduces to the Einstein-Maxwell scalar 
theory for 7 = 0. When 7=1 then (3.1) gives an action which is part of 
the low-energy action of string theory. For 7 = a/3 we get the action for 
the Kaluza-Klein theory. By varying the action (3.1) we obtain the following 




(3.1) 
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equations of motion: 



V, (e-2^*F*^) = 0, 

Rik = 2Vi^Vk^ + 2e-''^^Fi,F^-^gike-'^^F\ (3.2) 

Garfinkle, Horowitz and Strominger [27] obtained static spherically symmet- 
ric asymptotically flat black-hole solution described by the line element 

-(1 - —y-''r\de^ + sin'edct)'') (3.3) 
r 

with the dilaton field $ given by 

^(i-y)". (3.4) 



e 



and the electromagnetic field tensor component 



Ftr = % (3.5) 



where 

1-7^ 



1 + 7'' 

r_ and r+ are related to mass M and charge Q parameters as follows: 



(3.6) 



2 

_r_ 

:i+7^)' 



Q' - 7^^- (3-7) 



Virbhadra[72] proved that for Q = the GHS solution yields the Janis- 
Newman-Winicour solution [36] to the Einstein massless scalar equations. 
Virbhadra, Jhingan and Joshi[73] showed that the Janis-Newman-Winicour 
solution has a globally naked strong curvature singularity. 

For 7 = the solution yields the standard Reissner-Nordstrom of the 
Einstein-Maxwell theory, but for 7 7^ the solution is qualitatively different. 



47 



Certain qualitative features of the solutions of non-rotating charged dila- 
ton black-holes arc independent of the parameter 7. For instance, the sur- 
face r = r_|_ is an event horizon for all values of 7. A number of interesting 
properties of charged dilaton black holes critically depend on the dimension- 
less parameter 7 which controls the coupling between the dilaton and the 
Maxwell fields. The maximum charge, for a given mass, that can be carried 
by a charged dilaton black-hole depends on 7 [32]. When 7 7^ 0, the surface 
r = r_ is a curvature singularity while at 7 = the surface r = r_ is a non- 
singular inner horizon [33] . Both the entropy and temperature of these black 
holes depend on 7 [32]. The gyromagnetic ratio, i.e. the ratio of the magnetic 
dipole moment to the angular momentum, for charged slowly rotating dila- 
ton black holes depends on parameter 7 [32]. Chamorro and Virbhadra [12] 
showed, using Einstein's prescription, that the energy distribution of charged 
dilaton black holes depends on the value of 7. 



3.2 Energy distribution in charged dilaton black 
holes 

Virbhadra and Parikh[67] calculated, using the energy-momentum complex 
of Einstein, the energy distribution with stringy charged black holes (7 = 1) 
and found that 

Thus the entire energy of a charged black-hole in low-energy string theory is 
confined to the interior of the black-hole. Wc[77] found the same result using 
the Tolman definition. For the Reissner-Nordstrom metric, several definitions 
of energy give 

2r 

(Tod [59]; Hayward [31]; Aguirregabiria et al. [1]). Thus the energy is both 
in its interior and exterior. Chamorro and Virbhadra [12] studied, using the 
energy-momentum complex of Einstein, the energy distribution associated 
with static spherically symmetric charged dilaton black holes for an arbi- 
trary value of the coupling parameter 7 which controls the strength of the 
dilaton to the Maxwell field. We [78] investigated the energy distribution 
in the same spacetime in Tolman's prescription and got the same result as 
obtained by Chamorro and Virbhadra. The energy distribution of charged 
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dilaton black holes depends on the value of 7 and the total energy is in- 
dependent of this parameter. In this chapter we present the computations 
of the energy distribution for the Garfinkle-Horowitz-Strominger spacetime 
performed using the Tolman energy-momentum complex^. 



3.3 Tolman energy distribution 

We start by transforming the hne element (3.3) to quasi-Cartesian coordi- 
nates: 



r 1 



(1- 
<^-fl- 



-{xdx + ydy + zdz)"^, (3.8) 



according to 



9 





= COS 



^/x^~+y^~+^ 

--'(f)- 

Tolman's[61] energy-momentum complex is 



(3.9) 



TJ = —U, 



where 



^ kp 



9k^ap 



ak) 



(3.10) 



(3.11) 



^Virbhadra[74] pointed out that though the Tolman energy-momentum complex differs 
in form from the energy-momentum complex obtained by Einstein, both are equivalent in 
import. Therefore, the Tolman energy-momentum complex should be correctly referred to 
as the Tolman form of Einstein's energy-momentum complex. The author was not aware 
of this fact at the time of writing papers [77, 78] . 
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is the energy density, 7^"^ are the components of energy current density, 7^ 
are the momentum density components. Therefore, energy E for a stationary 
metric is given by the expression 

E=^J J J t/o% dxdydz. (3.12) 

After applying the Gauss theorem one has 

^ = ^ / / K"l^adS, (3.13) 

where /ia = {x/r,y /r, z/r) are the three components of a normal vector over 
an infinitesimal surface element dS = r^sin6d6d(j). 
The determinant of the metric tensor is given by 

«--(—) (3.14) 
and its non-vanishing contravariant components are: 

„<T+1 



r 



(r — (r 



- 



11 \^ ^ — / r 2 / \ AT 
9 — + ~ ^+^- J ~ ^ J ) 

12 ('■ - i-T'^ \ t ^ \ 1 

^ ^ ^:^T3 [r (r+ + r_) -r+r_J , 

22 ~ ''"-) [2/1 \ 41 

9 = — ;:^T3 — ly *^^^+ + ~ ^+^-'' - ^ J ) 

( \<J~l 

QQ ^ — / r 2 / \ 41 

g — \_z (rr+ + rr_ — r+r_j — r J , 

To compute the energy distribution using Eq. (3.13) we also require the 
following list of nonvanishing components of the Christoffel symbol of the 
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second kind. 



= x{Ci+C2X^), 

= y(ci + c2i/^), 

Tss = z{ci + C2z'^), 

= a;(c3 + C2|/^), 

^33 = a;(c3 + C2^^), 

r?! = ;2(a3 + a2a;^), 

= Z/(C3 + C2^^), 

= ^(C3 + C2y'), 

rl2 = y(c4 + C2a:;^), 

^13 = 2;(C4 + C2X^), 

r^i = a:;(c4 + C2y^), 

r23 = z{c4 + C2y'^), 

= X(C4 + C22;^), 

= y(c4 + c22;^), 

^00 ~ -^"^s ) 

T-i2 

Toi = XC6, 

^03 — ^^6, 

rL = = r?2 = C2xyz, (3.16) 

where 
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1 



Cl 



C2 = 



-[2r+r^ + 3r_r^ — 3r_r+r 



2H(r - r_) ' 
+r_^r+ — r_^r + (r_r+ — r_r — r+r — r^)r_(7], 
2r+^r + 6r_r+r — 3r_r^ — 3r+r^ — 3r_r+^ 



1 
2^ 



+ 



r — r+ 

(2r_ — r_(T)(r_r+ — r_r — r+r) 



C3 
C4 

C6 



2r4 
1 

2^ 



[r_r + 2r+r — r_r+ + (r — r+)r_(T], 



(r — r_) 



r — r_ 

2a-l 



2^2o-+4 



(r — r+)[(r — r_|_)r+ + (r — r+)r_(j], 



1 
2^ 



r+ r_(j 



(3.17) 



Using Eqs. (3.11) and (3.16) we obtain required components of Ul^ . These 



are 



X 



Uq^ — —T [r((Tr_ + r+) — (7r_r+] , 



(3.18) 



^0 



03 



—T [r((7r_ + r+) — (Tr_r+] . 



Now using (3.18) with (3.6) in (3.13) we get 



^(r) = M-^(l-f). 



(3.19) 



Thus, we get the same result as Chamorro and Virbhadra[12] obtained us- 
ing the Einstein energy-momentum complex. This is against the "folklore" 
that different energy-momentum complexes could give different and hence 
unacceptable energy distribution in a given spacetime. For the Reissner- 
Nordstrom metric (7 = 0) one gets E — M — (5^/2r, which is the same as 
obtained by using several other energy- momentum complexes (Aguirregabiria, 
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Figure 3.1: E /M on Z-axis is plotted against r/M on X-axis and 7 on Y-axis 
for Q/M = 0.1. 

Chamorro and Virbhadra) [1] and definitions of Penrose as well as Hayward 
[59], [31]. -E(r), given by (3.19), can be interpreted as the "effective gravita- 
tional mass" that a neutral test particle "feels" in the GHS spacetime. The 
"effective gravitational mass" becomes negative at radial distances less than 
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Chapter 4 

Energy Distribution in Ernst 
Space-time 

4.1 Introduction 

The well-known Melvin's magnetic universe[42] is a solution of the Einstein- 
Maxwell equations corresponding to a collection of parallel magnetic lines 
of force held together by mutual gravitation. Thorne[58] studied extensively 
the physical structure of the magnetic universe and investigated its dynam- 
ical behaviour under arbitrarily large radial perturbations. He showed that 
no radial perturbation can cause the magnetic field to undergo gravitational 
collapse to a space-time singularity or electromagnetic explosion to infinite 
dispersion. We[81] investigated the energy distribution in Melvin's magnetic 
universe and found encouraging results. The energy-momentum complexes 
of Einstein, Landau and Lifshitz, and Papapctrou give the same and accept- 
able energy distribution in Melvin's magnetic universe. A discussion of the 
Melvin's magnetic universe together with its energy distribution is given in 
sections (4.2) and (4.3), below. 

Ernst [24] obtained the axially symmetric exact solution to the Einstein- 
Maxwell equations representing the Schwarzschild black hole immersed in 
the Melvin's uniform magnetic universe. Virbhadra and Prasanna [69] stud- 
ied the spin dynamics of charged particles in the Ernst space-time. We [79] 
investigated energy distribution in the Ernst space-time and calculated the 
energy distribution using the Einstein energy-momentum complex. The first 
term of the energy expression is the rest-mass energy of the Schwarzschild 
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black hole, the second term is the classical value for the energy of the uniform 
magnetic field and the remaining terms in the expression are due to the gen- 
eral relativistic effect. The presence of the magnetic field is found to increase 
the energy of the system. Both the Ernst solution and a discuscussion of 
energy associated with a Schwarzschild black hole are given in sections (4.4) 
and (4.5). 



4.2 Melvin's magnetic universe 

The Einstein-Maxwell equations are 

where the energy-momentum tensor of the electromagnetic field is 



47r 



^ im-^ ' ^ ill ^ mn-^ 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



i^j*^ is the Ricci tensor and J* is the electric current density vector. 

Melvin [42] obtained the electrovac solution (J* — 0) to these equations 
which is expressed by the line element 



ds^ = [dt^ - dr^ - rV] - A-V^ sin^ 
and the Cartan components of the magnetic field are 

Hr = Bo COS 9, 



(4.5) 



where 



He = -A-^BoSinO, 

A = 1 + ^By sin^ e. 

4 ° 



(4.6) 



(4.7) 



Bo (= BovG/c^) is the magnetic field parameter and this is a constant in 
the solution given above. 
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The above solution had been obtained earher by Misra and Radhakrishna^ 
[44]. This space-time is invariant under rotation about, and translation along, 
an axis of symmetry. This is also invariant under reflection in planes com- 
prising that axis or perpendicular to it. Wheeler[76] demonstrated that a 
magnetic universe could also be obtained in Newton's theory of gravitation 
and showed that it is unstable according to elementary Newtonian analysis. 
However, Melvin[43] showed his universe to be stable against small radial 
perturbations and Thorne [58] proved the stability of the magnetic universe 
against arbitrary large perturbations. Thorne [58] further pointed out that 
the Melvin magnetic universe might be of great value in understanding the 
nature of extragalactic sources of radio waves and thus the Melvin solution to 
the Einstein-Maxwell equations is of immense astrophysical interest. We [79], 
[81] computed energy distribution in Mclvin's universe using the dcflnitions 
of Einstein, Landau and Lifshitz, and Papapetrou. For this space-time we 
found that these definitions of energy give the same and convincing results. 
The energy distribution obtained here is the same for all these energy mo- 
mentum complexes. In the next section we give the computations of energy 
distribution for this spacetime performed using energy-momentum complexes 
of Landau and Lifshitz, and Papapetrou. 



4.3 Energy distribution in Melvin magnetic 
universe 

The non-zero components of the energy-momentum tensor are 

Bl{l-2 sin^ e) 



87rA4 

Bl 



rp rp 3 

° ~ ' ~ 87rA4' 
^1 ^9 sin 9 cos 9 , , , 

= = ■ 

To get meaningful results using these energy-momentum complexes one is 
compelled to use "Cartesian" coordinates (see [46] and [74]). We use the 

^Melvin [43] mentioned that this solution is contained imphcitly as a special case among 
the solutions obtained by Misra and Radhakrishna. 
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following transformation: 



e = cos-^ 



= tan-i(^). (4.9) 
The line element (4.5) in t,x,y,z coordinates becomes 

ds' = A'df - A\dx' + dy' + dz') +(a^ + 1-] i^dy -yd^xY ^ .^^^ 



\ y x"^ + y"^ 

The determinant of the metric tensor is given by 

9 = -A", (4.11) 
and the non-zero contravariant components of the metric tensor are 



9'' 



5 

A-^x^ + A^y"^ 



x'^ + y"^ 

9'' = -(A 



2 ^ \ xy 



g22 



A-^y2 + A^x^ 



x^ + 

^33 ^ _^-2 (4 22) 

4.3.1 The Landau and Lifshitz energy- momentum com- 
plex 

The symmetric energy-momentum complex of Landau and Lifshitz [40] may 
be written as 

where 

^^'^■' = -^(^V (4.14) 
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is the energy density and are the momentum (energy current) density 
components, f'"-'"'^' has symmetries of the Riemann curvature tensor. The 
energy E is given by the expression 



^/ ji'^-'inpdS, (4.15) 



where /ip is the outward unit normal vector over an infinitesimal surface 
element dS. In order to calculate the energy component for Melvin's uni- 
verse expressed by the line element (4.10) we need the following non-zero 
components of t^^'' 



£0101 

^0102 
0202 



xy{K'^ - 1) 

^,2 _j_ ^2 ' 

1/2 + x^A*^ 



9 I 9 ' 

x-^ -\- 

= -1. (4.16) 
Equation (4.13) with equations (4.14) and (4.16) gives 

^00 ^ J_B^j^\ (4.17) 

For a surface given by parametric equations a; = r sin ^ cos 0, y = r sin sin 0, 
z = rcosO (where r is constant) one has ^/s = {x/r, y/r, z/r} and dS = 
r'^sin9d9d4>. Using equations (4.16) in (4.15) over a surface r = constant, we 
obtain 

Ell = -By + -By + —By + -^By. (4.18) 
6 " 20 ° 140 " 2520 " ^ ' 

4.3.2 The Energy-momentum complex of Papapetrou 

The Papapetrou energy-momentum complex[51]: 

(4.19) 

where 

M^^i^i = {g^'v''^ - ^'S'^ + g'^W' - gV) , (4.20) 
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is also symmetric in its indices, as discussed in chapter 2. are the energy 
and momentum density components. The energy E for a stationary metric 
is given by the expression 

^- = ^1 /-^""'-/^ ^^d^- (4.21) 

To find the energy component of the fine element (4.10), we require the 
following non-zero components of jV*-''^' 



»rooi2 ^ a;y(A^ - 1) 

9 I 9 ' 

jV°°33 ^ _2. (4.22) 
Equations (4.22) in equation (4.19) give the energy density component 

qOO ^ ±^2^3 (4 23) 

Thus we find the same energy density as we obtained in Section 4.3.1. We 
now use Eq. (4.22) in (4.21) over a 2-surface (as in the last Section) and 
obtain 

Ep = -By + -By + —By + -^By. (4.24) 

6 " 20 " 140 " 2520 ° ^ ' 

This result is expressed in geometrized units {G = 1 and c = 1). In the 
following we restore G and c and get 



Ep = -By + --,By + —^By + J—^lsy. (4.25) 

6 " 20 ° 140 c8 ° 2520 c^^ ° ^ ^ 



The first term — |— is the known classical value of energy and the rest of 
the terms are general relativistic corrections. The general relativistic terms 
increase the value of energy. 
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4.3.3 Discussion of Results 



In this section we obtained the energy distribution in Melvin's magnetic uni- 
verse. We used the energy-momentum complexes of Landau and Lifshitz, and 
Papapetrou. Both definitions give the same results {L^^ = Q^^,Eil = Ep). 
The first term in the energy expression (see equations (4.18) and (4.24)) is the 
well-known classical value for the energy of the uniform magnetic field and 
the other terms are general relativistic corrections. The general relativistic 
corrections increase the value of the energy. 



4.4 The Ernst solution 

Ernst [24] obtained an axially symmetric electrovac solution to the Einstein- 
Maxwell equations (4.1) to (4.4), describing the Schwarzschild black hole in 
Melvin's magnetic universe. The space-time is 



1 )df - (1 y^dr^ - r^de^ 

r r 



A-Vsin^ ed(f\ 



(4.26) 

with A given by (4.7) and the only Cartan component of the magnetic field 
which differs from those given in (4.6) is 

He = -A^'^Bo (1 - 2M/r)^/^ sin 6. (4.27) 

M and Bo arc constants in the solution. The non-zero components of the 
energy-momentum tensor are 



T ^ — —T 



2 _ Bl (2M sin^ ^ - 2r sin^ O + r) 



_ Bl{2Msm'e-r) 

^2 ^1 252(2M-r)sin^cos^ 

= -7^ = — — A ■ (4.28 

The Ernst solution is a black hole solution (r = 2M is the event horizon). 
For So = it gives the Schwarzschild solution and for M = it gives the 
Melvin's magnetic universe. The magnetic field has a constant value Bg 
everywhere along the axis. Ernst pointed out an interesting feature of this 
solution. Within the region 2M « r « B^^, the space is approximately 
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flat and the magnetic field approximately uniform, when |SoM| << 1. If the 
magnetic field is strong, i.e. |-Bo^| is of the order unity, then it tends to be 
more concentrated near the poles 9 = and 9 = 7:. 



4.5 The energy associated with Schwarzschild 
black hole in a magnetic universe. 

In this section we obtain the energy distribution of this spacetime using Ein- 
stein's energy-momentum complex. As in the case of the energy-momentum 
complexes of Landau and Lifshitz, and Papapetrou, to get meaningful re- 
sults for energy distribution in the prescription of Einstein one is compelled 
to use "Cartesian" coordinates. The line element (4.26) is easily transformed 
to "Cartesian" coordinates t,x,y,z using the standard transformation (4.9). 
We get 



A^{1 )dt^ 



- A2 
- 



A^ 



ax 



X 



+ A-2 



_|_ y2 



dx 



1 + 



r2(r - 2M) 

4:MXZ 



dz"- 



^2 |^2ax|/ ^ , ,_2 



r2(r - 2M) 



dxdz — A 



,2 

AMyz 



+ A 

dydz, 



2xy 



x^ + 



where 



2M 



a = 



r — 2M x^ + y 
Using the Einstein energy-momentum complex 

1 



2 ■ 



where 



/if = -K- 



Ik 



9ir. 



167r 



,ki 



'-g 



[-9 (/V™-^V™)]. 



the energy E for a stationary metric is given by the expression 

E=^J J J hl^dxdydz, 



dxdy 
(4.29) 

(4.30) 

(4.31) 
(4.32) 

(4.33) 
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and after applying the Gauss theorem, one has 

^ = / / h'oyads. 

4.5.1 Calculations 

The determinant of the metric tensor is given by 

9 = -A^ 

The non-zero contravariant components of the metric tensor are 



9 



00 _ 

9'' = 

9'' = A-^ 



r-2M' 



X 



2Mxy 



xy 



x"^ + y"^ 



-A' 
+ A2 



x'^ + y'^ 



xy 



x"^ + y"^ 



9 



22 



-A" 



9'' = A-2 
/3 = A-2 



2M7/2 
2Mz2 

2Mxz' 
2Myz 



x^ + 1/^ 



-A2 









(4.34) 



(4.35) 



(4.36) 



The only required components of h^^ in the calculation of energy are the 
following: 



u03 



4Mx 



+ (A^ - 1) 



X 



x"^ + y"^ 

y 

x^ + 1/^ 



(4.37) 
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Now using (4.37) with (4.34) we obtain the energy distribution in the Ernst 
space-time. 

1 i*'2'K 

E^M+ — / (A^-l)rsined^#. (4.38) 

16vr Jg^o .U=o 

We substitute the value of A in the above and then integrate. We get 

E^M + -By + —By + —By + -^By. u.m 

6 ° 20 " 140 " 2520 ° ^ ' 

The above result is expressed in geometrized units (gravitational constant 
G = 1 and the speed of light in vacuum c — 1). In the following we restore 
G and c and get 

E . Mc^ . \By . l^py ^ ±^py + ^^l^ny. (4.4o) 

The first term Mc^ is the rest-mass energy of the Schwarzschild black hole, 
the second term ^By is the well-known classical value of the energy of the 
magnetic field under consideration, and the rest of the terms are general rel- 
ativistic corrections. For very large B^r, the general relativistic contribution 
dominates over the classical value for the magnetic field energy. As men- 
tioned in Section 2, the gravitational field is weak for 2M « r « B~^ (in 
G = 1, c = 1 units). Thus in the weak gravitational field we have Bqt << 1 
; therefore, the classical value for the magnetic field energy will be greater 
than the general relativistic correction in these cases. 

4.5.2 Discussion of Results 

In the above section we considered the Ernst space-time and calculated the 
energy distribution using the Einstein energy-momentum complex. It beauti- 
fully yields the expected result: The first term is the Schwarzschild rest-mass 
energy, the second term is the classical value for energy due to the uniform 
magnetic field {E = J J J B'^dV, where dV is the infinitesimal volume el- 
ement, yields exactly the same value as the second term of (4.40) ), and the 
rest of the terms are general relativistic corrections. The general relativistic 
terms increase the value of the energy. 
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4.6 Conclusion 



In this chapter we calculated energy distributions in Melvin's magnetic uni- 
verse and Ernst space-time using prescriptions of Landau and Lifshitz, Pa- 
papetrou, and Einstein. We got encouraging results for the asymptotically 
non-fiat space-times considered above. We also note that the results ob- 
tained indicate that the energy-momentum complexes of Einstein, Landau 
and Lifshitz, and Papapetrou give the same energy distribution for the Melvin 
magnetic universe. It was believed that the results are meaningful for en- 
ergy distribution in the prescription of Einstein only when the space-time 
studied is asymptotically Minkowskian. However, recent investigations of 
Rosen and Virbhadra [56], Virbhadra [70], Aguirregabiria et al. [1], and 
Xulu [79], [81] showed that many energy- momentum complexes can give 
the same and appealing results even for asymptotically non-flat space-times. 
Aguirregabiria et al. showed that many energy-momentum complexes give 
the same results for any Kerr-Schild class metric. There are many known 
solutions of the Kerr-Schild class which are asymptotically not fiat. For ex- 
ample, Schwarzschild metric with cosmological constant. The general energy 
expression for any Kerr-Schild class metric obtained by them immediately 
gives E — M — {X/3)r^ where A is the cosmological constant. This result 
is very much convincing. A > gives repulsive effect whereas A < gives 
attractive effect. 
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Chapter 5 



Total Energy of the Bianchi 
Type I Universes 

5.1 Introduction 

A wide range of cosmological models may be deduced from Einstein's field 
equations. The 1965 observation of the cosmic microwave background radia- 
tion, by Penzias and Wilson, was the most important cosmological discovery 
since Bubble's 1929 announcement that all galaxies recede from us at veloc- 
ities proportional to their distance from us. This discovery of background 
radiation strongly support that some version of the big bang theory is cor- 
rect and it also resulted in some conjectures regarding the total energy of 
the universe. Tryon[62], assuming that our Universe appeared from nowhere 
about 10^° years ago, remarked that the conventional laws of physics need 
not have been violated at the time of creation of the Universe. He proposed 
that our Universe must have a zero net value for all conserved quantities. 
The arguments he presented indicate that the net energy of our Universe 
may be indeed zero. His big bang model (in which our Universe is a fluctua- 
tion of the vacuum) predicted a homogeneous, isotropic and closed Universe 
consisting of matter and anti-matter equally. Tryon [62] also cited an elegant 
topological argument by Bcrgmann that any closed universe has zero energy. 

Two decades later, the work of Cooperstock [18] and Rosen [55] revived 
the interest in the investigations of the energy of the Universe. Cooperstock 
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[18] considered the conformal Priedmann-Robertson- Walker (FRW) metric 



dr"^ 



kr'^ 



(5.1) 



By making use of calculations involving killing vectors, Coopcrstock and 
Israelit (see in [18]) were able to express the covariant conservation laws 
T**^ ;fc — 0, in the form of an ordinary divergence: 



= 0. 



(5.2) 



^From (5.2), Cooperstock [18] was able to conclude that the total density of 
the universe is zero. 

Rosen [55] considered a closed homogeneous isotropic universe described 
by the Friedmann-Robertson- Walker (FRW) metric: 



ds^ = de 



a{ty 



{dr"^ + r'^de'^ + r'^sm^edip^) . (5.3) 



Then using Einstein's prescription, he obtained the following energy-momentum 
complex^ 



e°- — 



, . (5.4) 
.(l + rV4)' (l + rV4)^ ^ ' 

By integrating the above over all space, one finds that the total energy E 
of the imiverse is zero. These interesting results fascinated some general 
relativists, for instance, Johri et al. [37], Banerjee and Sen [3] and Xulu [80]. 
Johri et al. [37], using the Landau and Lifshitz energy-momentum complex, 
showed that the total energy of an FRW spatially closed universe is zero at all 
times irrespective of equations of state of the cosmic fluid. They also showed 
that the total energy enclosed within any finite volume of the spatially fiat 
FRW universe is zero at all times. In this chapter we investigate the total 
energy of the Bianchi type I universes. 

^To avoid any confusion we mention that we use the term energy- momentum complex 
for one which satisfies the local conservation laws and gives the contribution from the 
matter (including all non-gravitational fields) as well as the gravitational field. Rosen [55] 
used the term pseudo-tensor for this purpose. We reserve the term energy-momentum 
pseudotensor for the part of the energy-momentum complex which comes due to the grav- 
itational field only. 
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5.2 Bianchi type I space-times 



The Bianchi type I space-times arc expressed by the hne element 

ds^ = dt-" - e^'dx" - e^^dy'' - e^^dz\ (5.5) 

where I, m, n are functions of t alone. The nonvanishing components of the 
energy-momentum tensor Tj ( = ^G^^, where G^^ is the Einstein tensor) 
are 

Tq° = — ( lfn + rhn -|- n/ V 

7\ ^ = — ( + 71^ + rhn + 111 + 71), 

Stt 

T^^ = — (^fi^ + i^ + ni + n + i^ , 

= ^ P + m^ + /m + [ + m ) . (5.6) 

The dot over l,m,n stands for the derivative with respect to the coordinate 
t. The metric given by Eq. (5.5) reduces to the spatially flat Friedmann- 
Robertson- Walker metric in a special case. With l{t) = m{t) = n{t), defining 
R{t) — e'^*^ and transforming the line element (5.5) to t,x,y,z coordinates 
according to x = r sin ^ cos 0, y = r sin ^ sin 0, z — r cos 6 gives 

ds'' = df - [R {t)f {dr^ + {de^ + sin^ed(j)^) }, (5.7) 

which describes the well-known spatially flat Friedmann-Robertson- Walker 
space-time. 



5.3 Energy distribution in Bianchi type I space- 
times 

The Bianchi type / solutions, under a special case, reduce to the spatially fiat 
FRW solutions. Banerjee and Sen [3] studied the Bianchi type / solutions, us- 
ing the Einstein energy-momentum complex, and found that the total (mat- 
ter plus field) energy density is zero everywhere. As the spatially flat FRW 
solution is a special case of the Bianchi type I solutions, one observes that 
the energy-momentum complexes of Einstein and Landau and Lifshitz give 
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the same result for the spatially flat FRW solutions. Because there is a per- 
ception that different complexes could give different and hence meaningless 
results for a given metric, we [80] investigated energy-momentum distribution 
in Bianchi type I space-times using energy-momentum complexes of Landau 
and Lifshitz, Papapetrou, and Weinberg. The results of these investigations 
are presented in this section. 



5.3.1 The Landau and Lifshitz energy- momentum com- 
plex 

In order to calculate the energy and momentum density components of the 
line element (5.5) using the symmetric energy-momentum complex of Landau 
and Lifshitz [40]: 

where 

S^3kl^_g^gi3gkl_gikgil^^ ^^^^^^ 

the required nonvanishing components of S'^^^^ are 

^0101 _ _g2m+2n 

^0110 _ g2m+2n 

^0202 _ _g2i+2n 

^0220 _ g2Z+2n 

^0303 _ _g2«H-2m 

^0330 ^ g2Z+2m ^5 ^^q) 

Using the above results in (5.8) and (5.9) we obtain the energy density and 
energy current (momentum) density components, respectively as: 

= = 0. (5.11) 

Hence the energy and momentum components 

U^dx^dx^dx^ (5.12) 

vanish. 
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5.3.2 The Energy-momentum complex of Papapetrou 

In order to calculate the energy and momentum density components of the 
line element (5.5) using the symmetric energy-momentum complex of Papa- 
petrou [51]: 

^ il^^"'i (^-'^^ 

where 

M'^^^ = [g'^T]^^ - g'^r]^^ + g^^r]'^ - g^^rf^) (5.14) 

and rf^ is the Minkowski metric, we require the following nonvanishing com- 
ponents of jV'^'^': 





= -(1 + e- 


-2l^^l+m+n 




^—l+m+n 




jy0022 


= -(1 + e- 


-2m\ l+m+n 


jy^0220 


l—rn+n 




jy0033 


= -(1 + e- 


-2n^gi+m+n 


jy^0330 


l+m—n 




_y^0101 


_ ^0202 _ 





(5.15) 

Using the above results in (5.13) and (5.14) wc obtain the energy density 
component and momentum (energy current) density components^]"'' as: 

= = 0. (5.16) 

Again, we find that the energy and momentum components 

n'^'dx^dx^dx^ (5.17) 

J J J 
vanish. 

5.3.3 The Weinberg energy-momentum complex 

The energy and momentum density components of the line element (5.5) were 
calculated using the symmetric energy- momentum complex of Weinberg [75] : 

VF'^ = — A'^\ (5.18) 
167r ^ ^ 
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where 

- off - - 0^^' + + ^ - ^'-''^ 

and 

hij = g,j - rjij- (5.20) 

rjij is the Minkowski metric. In this case, using the equations (5.5) and 

(5.19), we find that all the components of A*-''^ vanish. Thus Eq. (5.18) 
yields 

W''' = 0. (5.21) 
Therefore the energy and momentum components 

W'^dx^dx'^dx^ (5.22) 

J J J 

also vanish. 




5.4 Conclusion 

In recent years some researchers showed interest in studying the energy con- 
tent of the universe in different models (see Cooperstock [18] , Rosen [55], 
Johri et al. [37], Banerjee and Sen [3] . Cooperstock [18] investigated energy 
density for conformal Friedmann-Robcrtson- Walker metric and by making 
use of calculations involving killing vectors he was able to deduce that the 
total energy density is equal to zero. Rosen [55] studied the total energy of 
a closed homogeneous isotropic universe described by the FRW metric using 
the Einstein energy-momentum complex, and found that to be zero. Using 
the Landau and Lifshitz prescription of energy and momentum Johri et al. 
[37] demonstrated that (a) the total energy of an FRW spatially closed uni- 
verse is zero at all times irrespective of equations of state of the cosmic fiuid 
and (b) the total energy enclosed within any finite volume of the spatially 
flat FRW universe is zero at all times. Banerjee and Sen [3] showed that 
the energy and momentum density components vanish in the Bianchi type I 
space-times (they used the energy-momentum complex of Einstein) . 

It is usually suspected that different energy-momentum complexes could 
give different results for a given geometry. Therefore, we[80] extended the 
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investigations of Banerjee and Sen with three more energy-momentum com- 
plexes (proposed by Landau and Lifshitz, Papapetrou, and Weinberg) and 
found the same results (see equations (5.11), (5.16) and (5.21)) as reported 
by them. Note that the energy density component of the energy-momentum 
tensor is not zero for the Bianchi type I solutions (see Eq. (5.6)); however, it 
is clear from equations (5.11), (5.16) and (5.21) that the total energy density 
(due to matter plus field, as given by the energy-momentum complexes) van- 
ishes everywhere . This is because the energy contributions from the matter 
and field inside an arbitrary two-surface in Bianchi type I space-times can- 
cel each other. These results illustrate the importance of energy-momentum 
complexes (as opposed to the perception against them that different com- 
plexes could give different and hence meaningless results for a given metric) 
and also supports the viewpoint of Tryon. 
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Chapter 6 



M0ller Energy for the 
Kerr-Newman Metric 

6.1 Introduction 

The investigations of Hawking, Israel, Carter, Robinson and others on the 
properties of black holes built-up to the proof of the so-called "No Hair" the- 
orem which shows that black holes are completely described by only three 
quantities namely mass M, charge e, and angular momentum a (see in Israel 
[35]). Hence, the stationary axially symmetric and asymptotically flat Kerr- 
Newman solution which is parameterized by mass M, charge e, and angular 
momentum a, is the most general black hole solution to the Einstein-Maxwell 
equations. This solution describes the exterior gravitational and electromag- 
netic fields of a charged rotating object. When e = 0, it describes the Kerr 
family of axially symmetric solutions that give the geometry of space-time 
surrounding rotating uncharged objects. When a = 0, it describes the spher- 
ically symmetric Reissner-Nordstrom solution of charged non-rotating black 
holes. For both a = and e = the solution reduces to the spherically 
symmetric Schwarzschild solution of the simplest type of black hole which 
is only characterized by the mass M . The Kerr-Newman solution, as the 
most general black hole solution, is therefore of vital importance in studying 
the geometry surrounding compact objects. In this chapter we investigate 
energy distribution in Kerr-Newman space-time. 

The energy distribution in the Kerr-Newman (KN) space-time was ear- 
lier computed by Cohen and de Felice [16] using Komar's prescription. Virb- 
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hadra ([63], [64]) showed that, up to the third order of the rotation parameter, 
the energy-momentum complexes of Einstein and Landau-Lifshitz give the 
same and reasonable energy distribution in the KN space-time when calcu- 
lations are carried out in Kerr-Schild Cartesian coordinates. Cooperstock 
and Richardson [21] extended the Virbhadra energy calculations up to the 
seventh order of the rotation parameter and found that these definitions give 
the same energy distribution for the KN metric. Aguirregabiria et al. [1] 
performed exact computations for the energy distribution in KN space-time 
in Kerr-Schild Cartesian coordinates. They showed that the energy distri- 
bution in the prescriptions of Einstein, Landau-Lifshitz, Papapetrou, and 
Weinberg (ELLPW) gave the same result. In a recent paper Lessner[41] in 
his analysis of M0ller's energy-momentum expression concludes that it is a 
powerful concept of energy and momentum in general relativity. We [82] eval- 
uated the energy distribution in KN field using the M0ller energy-momentum 
prescription. The results of our investigation [82] are given below. 

In a series of papers ([17], [19], [20]), Cooperstock has propounded a hy- 
pothesis which essentially states that the energy and momentum in a curved 
space-time are confined to the regions of non-vanishing energy-momentum 
tensor of the matter and all non-gravitational fields. It is of interest to 
investigate whether or not the Cooperstock hypothesis holds good. Our re- 
sults ([82], [83]) and the recent results of Bringley [7] support this hypothesis. 
In this chapter we use the Kerr-Newman space-time for testing the Coop- 
erstock hypothesis. We first give the KN metric, followed by M0ller energy 
distribution and a discussion of results. 



6.2 The Kerr-Newman metric 

The Kerr-Newman metric in Boyer-Lindquist coordinates {t, p, 9, (f)) is ex- 
pressed by the line element: 

ds^ ^^[dt-a sin^ ^d0]2 - -^[(p2 + a^) d<j) - adtf - -j-dp" - r^dO^ (6.1) 
'"o '"o ^ 

where A := p^ — 2Mp-\-e'^-\-a? and Tq := p^-\-a^ cos^ 9. M, e and a are respec- 
tively mass, electric charge and rotation parameters and the corresponding 
electromagnetic field tensor is: 

F — erQ'^Kp'^ ~ cos^ 9)dp A dv — 20^ p cos 9d9 A dv 

- a sin^ 9{p'^ - cos^ 9)dp Ad(j) + 2ap{p^ + a^) cos 9 sin 9d9 Ad^]. 
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(For details see in Carter [10].) The KN space-time has {p — constant} null 
hypersurfaces for g'''' — 0, which are given by 

p± = M ± VM2 - e2 - a2 . (6.3) 

There is a ring curvature singularity p = in the KN space-time. This space- 
time has an event horizon ai p — p^. It describes a black hole if and only if 

M2 > e2 + a\ 

The Boycr-Lindquist coordinates are singular at p = p±. Therefore, to 
remove this coordinate singularity t is replaced with a null coordinate and (j) 
with an 'untwisting' angular coordinate (p using the following transformation: 

dt — dv — ^—T — dp, 
A 

a, 

d4> = dip — '^^Py (6-4) 

and thus we express the KN metric in advanced Eddington-Finkelstein co- 
ordinates (Misner et. al. [45] refer to these as Kerr coordinates) (v, p, 6*, ip) 
as: 

ds^ = (l - + % \ dv^ - 2dv dp + ^ (2Mp - e^) dv dip - rldQ 



2 



-\-2a sin^ Odpdip — 



(p2 + a^) sin2 9 + '^-^l^a' sin^ 



dip\ (6.5) 



The energy-momentum complexes of Einstein, Landau- Lifshitz, Papapetrou 
and Weinberg are coordinate-dependent and require the use of quasi-Cartesian 
coordinates. Thus we transform the above to Kerr-Schild Cartesian coordi- 
nates (T, X, y, z) according to: 

T ^ v-p, 

X — sin 6 {p cos (/? -|- a sin (/?) , 

y — sin ^ (p sin (/? — a cos </?) , 

z — pcosO, (6.6) 
and one has the line element 



ds^ = dT^ - dx^ - dy^ - dz" - ''^f X 

^ p4 + a?z^ 

2 



dT H — — - ixdx + ydy) H — — - (ydx — xdy) H — dz ] (6.7) 
-\- p^ + p^ p 
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The components of the energy-momentum tensor T^, in quasi-Cartesian co- 
ordinates, are given by 





+ 


2ay 


—2ax 







-2ay 




2xy 


2xz 


87rr6 


2ax 


2xy 


-{R^ + a^-2y^) 


2yz 





2xz 


2yz 





(6.8) 

where 

ro' = (i?2-a2)2 + 4aV (6.9) 
(for details see in Cooperstock and Richardson [21]). 



6.3 Energy distribution in Kerr-Newman met- 
ric. 

In this Section we first give the energy distribution in the KN space-time ob- 
tained by some authors and then using the M0ller energy-momentum complex 
we obtain the energy distribution for the same space-time. 



6.3.1 Previous results 



The energy distribution in Komar's prescription obtained by Cohen and de 
Felice[16], using the KN metric (6.1) in Boyer-Lindquist coordinates, is given 

by 



K 



M 

2p 



1 + 



ap 



arctan ( - 
P 



(6.10) 



(The subscript K on the left hand side of the equation refers to Komar.) 
Aguirregabiria et al\l] studied the energy- momentum complexes of Einstein, 
Landau-Lifshitz, Papapetrou and Weinberg for the KN metric. They showed 
that these definitions give the same results for the energy and energy current 
densities. They used the KN metric (6.7) in Kerr-Schild Cartesian coordi- 
nates. They found that these definitions give the same result for the energy 
distributon for the KN metric, which is expressed as 



-Eellpw — ^ ~ it 
Ap 



1 + ^ — arctan - 

op \p 



(6.11) 



75 



(The subscript ELLPW on the left hand side of the above equation refers 
to the Einstein, Landau-Lifshitz, Papapetrou and Weinberg prescriptions.) 
It is obvious that the Komar definition gives a different resuh for the Kerr- 
Newman metric as compared to those obtained using energy-momentum com- 
plexes of ELLPW. However, for the Kerr metric (e = 0) all these definitions 
yield the same results. These results obviously support the Cooperstock 
hypothesis. 



6.3.2 The M0ller energy distribution 

In order to calculate the energy and momentum density components of the 
Kerr-Newman metric using the M0ller energy- momentum complex [46] SJ^'^ 
given by 

%'-l^Xi'\i , (6.12) 
where the antisymmetric superpotential Xi^'' is 

Xi"' = -Xi'' - hn,m - 9im,n] (6.13) 

the only required non- vanishing component of Xi'^'' is 

Xo'' - yy {Ma' cos^ e-Mp' + e'p) . (6.14) 



1 ' ' 0/3 



Using the above expression in 

^^i^J J Xo' PpdS, (6.15) 

for the energy £^ of a stationary metric, (where //^ is the outward unit normal 
vector over an infinitesimal surface element dS), we then obtain the energy 
E inside a surface with {p — constant} given 



,2 



Em01 ^M- 

2p 



^ (a^ + p ) I a 



1 H arctan 

ap \P 



(6.16) 



(The subscript M0I on the left hand side of this equation refers to M0ller's 
prescription.) 
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6.3.3 Discussion of Results 



The above result (Eq. 7.5) , obtained using M0ller's complex, agrees with 
the energy distribution (Eq. 6.10) obtained by Cohen and de Felice[16] in 
Komar's prescription. It differs by a factor of two in the second term of 
the energy distribution from that (Eq. 6.11) computed by Aguirregabiria et 
al. using ELLPW complexes. However, in both cases the energy is shared 
by both the interior and exterior of the KN black hole. It is clear that 
the definitions of ELLPW, Komar, and now that of M0ller also upholds the 
Cooperstock hypothesis for the KN metric. The total energy (p oo in all 
these energy expressions) give the same result M. 




25 50 75 100 



X 



Figure 6.1: £^ellpw and 8km on Z-axis are plotted against IZ on X-axis and 
S on Y-axis for Q = 0.1. The upper (transparent one) and lower surfaces 
are for £^ellpw and Skm respectively. 

Now defining 

r> _ -^ELLPW c _ -^'ELLPW c _ "^K Em01 
i^ELLPW '■— 77 ; i^ELLPW •— 77 ; ^KM 



e p 



5:=^, Q:=^, 7^:=^ (6.17) 
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the equations (6.10), (6.11) and (7.5) may be expressed as 



and 



^■KM — 1 



27^ 



1 + 



s n\ fs 

+ ^ arctan — 

n s \n 



(6.18) 



(6.19) 



The ring curvature singularity in the KN metric is covered by the event 
horizon for (Q^ + 5^) < 1 and is naked for (Q^ + «S^) > 1. In Fig. 1 we plot 
^^ELLPW and £^KM against IZ and iS for Q = 0.1. As the value of TZ increases 
the two surfaces shown in the figure come closer. 



6.4 Conclusions 

The prescriptions of Einstein, Landau and Lifshitz, Papapctrou, and Wein- 
berg used for calculating the energy-momentum distribution in a general 
relativistic system restrict one to make calculations in quasi-Cartesian co- 
ordinates. This shortcoming of singling out a particular coordinate system 
prompted M0ller[46] to construct an expression which enables one to eval- 
uate energy in any coordinate system. According to M0ller[46] this expres- 
sion should give the same values for the total energy and momentum as 
the Einstein's energy-momentum complex for a closed system. However, 
M0ller's energy-momentum complex was subjected to some criticism (see in 
M0ller[47], Kovacs[39], Novotny[50]). Further Komar[38] formulated a new 
definition of energy in a curved space-time. This prescription, though not 
restricted to the use of "Cartesian coordinates", is applicable only to the 
stationary space-times. The M0ller energy-momentum complex is neither 
restricted to the use of particular coordinates nor to the stationary space- 
times. Recently, Lessner[41] pointed out that the M0ller definition is a pow- 
erful concept of energy and momentum in general relativity. However, it is 
worth noting that for the Reissner- Nordstrom metric -Bellpw — M — e^/ (2p) 
(the Penrose definition also gives the same result [74] and this provides the 
weak field limit) whereas i?KM0i = M — / p does not give the weak field 
limit. This question must be investigated carefully. In the next chapter we 
investigate the Cooperstock hypothesis for the non-static space-times with 
M0ller's energy-momentum complex. 
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Chapter 7 



Energy of the Nonstatic 
Spherically Symmetric Metrics 

7.1 Introduction 

In a recent paper, Virbhadra[74] investigated whether or not the energy- 
momentum complexes of Einstein, Landau and Lifshitz (LL) , Papapetrou, 
and Weinberg give the same energy distribution for the most general non- 
static spherically symmetric metric and, contrary to previous results of many 
asymptotically flat spacetimes [1, 12, 13, 21, 63, 64, 65, 66, 67, 68, 70, 72, 
77, 78] and asymptotically non-flat spacetimes [56, 70, 79, 80], he found that 
these definitions disagree. He observed that the energy-momentum complex 
of Einstein gave a consistent result for the Schwarzschild metric whether one 
calculates in Kerr-Schild Cartesian coordinates or Schwarzschild Cartesian 
coordinates. The prescriptions of LL, Papapetrou and Weinberg furnish the 
same result as in the Einstein prescription if computations are carried out 
in Kerr-Schild Cartesian coordinates; however, they disagree with the Ein- 
stein definition if computations are done in Schwarzschild Cartesian coordi- 
nates. Thus, the definitions of LL, Papapetrou and Weinberg do not furnish 
a consistent result. Based on this and some other investigations (see also in 
Bergqvist [5] , Bernstein and Tod [4], Virbhadra concluded that the Einstein 
method seems to be the best among all known (including quasi-local mass 
definitions) for energy distribution in a space-time. 

In the previous chapter we highlighted Lessner's arguments indicating the 
importance of the M0ller energy-momentum expression. So in the present 



79 



chapter we wish to revisit the M0ller energy- momentum prescription by 
presenting the resuh of our investigation [83] of the energy distribution in 
the most general nonstatic spherically symmetric space-time using M0ller's 
energy-momentum complex. This result is compared with the Virbhadra 
energy expression obtained by using the energy-momentum complex of Ein- 
stein. We also discuss some examples of energy distributions in different 
prescriptions. In the next section we give the energy expression obtained by 
Virbhadra [74]. 



7.2 Virbhadra's energy result 

The most general nonstatic spherically symmetric space-time is described by 
the line element 

ds"^ = a{r, t) dt^-(5{r, t) dr^-2-f{r, t) dt dr-a{r, t) {dO'^ +sm^ dcp^) . (7.1) 

This has, amongst others, the following well-known space-times as special 
cases: The Schwarzschild metric, Reissner-Nordstrom metric, Vaidya met- 
ric, Janis-Newman-Winicour metric, Garfinkle-Horowitz-Strominger metric, 
a general non-static spherically symmetric metric of the Kerr-Schild class 
(discussed in Virbhadra's paper[74]). Virbhadra[74] explored the energy dis- 
tribution in the most general nonstatic spherically symmetric space-time 
(7.1) using the energy-momentum complex of Einstein (2.16). To com- 
pute the energy E (= Pq) using the Qq^ component of the Einstein energy- 
momentum complex Virbhadra transformed the line element (7.1) to "Carte- 
sian coordinates" (t, a;, z) using a; = r sin ^ cos (p^y = rsm.6 sin (f),z = r cos 9 
and t remaining the same. Then using (2.20) he obtained the energy distri- 
bution which is given below: 

P r[a{P-cr- ra ^) - 7 (ra ^ - 7)] 

-C'Einst — / ^ [ I 

where comma indicates partial differentiation. In the next Section we obtain 
the energy distribution for the same metric in M0ller's formulation. 
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7.3 Energy distribution in M0ller's formula- 
tion 



In this Section we use the energy-momentum complex of M0ller to obtain 
energy distribution in the most general nonstatic spherically symmetric met- 
ric given by the equation (7.1). Since the M0ller complex is not restricted to 
the use of "Cartesian coordinates" we perform the computations in t, r, 9, (f) 
coordinates, because computations in these coordinates are easier compared 
to those in t, x, y, z coordinates. 

The M0ller energy- momentum complex '^■'^ is given by (2.47), with the 
anti-symmetric superpotential Xj*^' given by (2.48). To compute the energy 
distribution 

^ = ^ / / Xo^l^pdS, (7.3) 
for the line element (7.1) under consideration we calculate 

XO — , ,1/2 ' \'-^) 

which is the only required component of Xi^ for our purpose. 

Using the above expression in equation (6.15) we obtain the energy dis- 
tribution 

_ (a,r + l,t) (rr^ 
2 [afi + 7^) 

It is evident that the energy distribution for the most general nonstatic spher- 
ically symmetric metric the definitions of Einstein and M0ller disagree in gen- 
eral (compare (7.2) with (7.5)). However, these furnish the same results for 
some space-times, for instance, the Schwarzschild and Vaidya space-times [66]. 
In the next Section we will compute energy distribution in a few space-times 
using (7.2) and (7.5). 



7.4 Examples 



In this Section we discuss a few examples of space-times in the Einstein as 
well as the M0ller prescriptions. We also test the Cooperstock hypothesis 
with these examples. 
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1. The Schwarzschild solution 

This solution is expressed by the line element 

ds^ =(^1- dt' - (^1 - ' dr^ - {de + sin^ eHct>'') . 

(7.6) 

Equations (7.2) and (7.5) furnish (see also in [46, 72]) 

E-Einst — Em0\ — M (7-7) 

showing that these two definitions of energy distribution agree for the 
Schwarzschild space-time and the above results support the Cooper- 
stock hypothesis. 

2. The Reissner-Nordstrom solution 

The Reissner-Nordstrom solution is given by 

ds^^ — + de- — + ' dr'-r' (d0' + sin' ed<f>') , 

(7.8) 

and the antisymmetric electromagnetic field tensor 

F,. = ^ , (7.9) 

where M and e are respectively the mass and electric charge parame- 
ters. 

For this space-time equations (7.2) and (7.5) furnish (see also in [64, 21]) 

^^Einst^M-^ (7.10) 

zr 

and 

Em01 = M--. (7.11) 

r 

Both of these results obviously support the Cooperstock hypothesis. 

3. The Janis-Newman-Winicour solution 

This solution has been usually incorrectly referred to in the literature 
as the Wyman solution. Virbhadra[72] proved that the Wyman solu- 
tion is the same as the Janis-Newman-Winicour solution. As Janis, 
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Newman and Winicour obtained this solution much before Wyman, 
Virbhadra[74] rightly referred to this as the Janis-Newman- Winicour 
solution. This solution is given by 

ds^ ^{l-^y dt'- (l - f ) ' dr'- \^ {d9^ + sin^ 9 d^') 

(7.12) 

and the scalar field 

^^4^1n fl- (7.13) 

where 



2M 

^ = -B- 

B = 2 + g2. (7.14) 

M and q are the mass and scalar charge parameters respectively. For 
q — this solution furnishes the Schwarzschild solution. 

Virbhadra[74] computed the energy expression for this metric using Eq. 
(7.2). We do the same here using equation (7.5). Thus we find that 

^Einst = ^M0l = M, (7.15) 

which shows that these two definitions of energy distribution agree for 
the Janis-Newman- Winicour space-time. 

4. Garfinkle-Horowitz-Strominger solution 

The Garfinkle-Horowitz-Strominger static spherically symmetric asymp- 
totically fiat solution (see in [27]) is described by the line element (3.3). 
In order to compute the energy distribution in Garfinkle-Horowitz- 
Strominger space-time using the energy-momentum complex of Ein- 
stein, Chamorro and Virbhadra[12] transformed (3.3) to quasi-Cartesian 
coordinates (3.8) . Then, by making use of (2.16) they found the fol- 
lowing expression in Einstein prescription: 

£^Einst = M - 1^ (1 - A^) . (7.16) 
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We compute the energy distribution for the Garfinkle-Horowitz-Strominger 
space-time using equation (7.5) and obtain 



Em01^M--{1-X^) . (7.17) 

Thus, these two definitions give difi^erent results (there is a difference 
of a factor 2 in the second term) but they obviously support the Coop- 
erstock hypothesis. 

Now defining 

the equations (7.16) and (7.17) may be expressed as 

W = 1-^(1-A2) (7.19) 

and 

Sm01 = 1 - (1 - (7-20) 

For = 1, £^Einst = ^M0i = M; however, they differ for any other values 
of A^. For any values of A^ < 1, ^Einst as well as £m0\ decrease with 
an increase in and increase with increase in TZ. SEmst > ^M0\ and 
they asymptotically {TZ oo) reach the value 1. The situation is just 
opposite for any values of A^ > 1 : £^Emst as well as Sm0\ increase with 
an increase in and decrease with increase in R. ^Einst < ^M0i and 
they asymptotically {TZ — > oo) reach the value 1. 

We plot the energy distributions £^Einst and Sm01 for A = (Reissner- 
Nordstrom space-time) in the figure 1 and for A^ = 1.2 in figure 2. 
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0.1 0.075 0.05 0.025 

Y 

Figure 7.1: £^Einst and Suoi on Z-axis are plotted against 71 on X-axis and Q 
on Y-axis for A = (Reissner- Nordstrom metric). The upper (grid-like) and 
lower surfaces are for £^Einst and Sm01 respectively. 

7.5 Conclusion 

Based on some analysis of the results known with many prescriptions for en- 
ergy distribution (including some well-known quasi-local mass definitions) in 
a given space-time Virbhadra[74] remarked that the formulation by Einstein 
is still the best one. Lessner[41] argued that the M0ller energy-momentum 
expression is a powerful concept of energy and momentum in general rel- 
ativity, which motivated us to study this further. We obtained the energy 
distribution for the most general nonstatic spherically symmetric metric using 
M0ller's definition. The result we found differs in general from that obtained 
using the Einstein energy-momentum complex. However, these agree for the 
Schwarzschild, Vaidya and Janis-Newman-Winicour space-times. They dis- 
agree for the Reissner-Nordstrom space-time. For the Reissner-Nordstrom 
space-time ^'Einst = M — e^/ (2r) (the seminal Penrose quasi-local mass def- 
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Figure 7.2: £^Einst and £m01 on Z-axis are plotted against 71 on X-axis and Q 
on Y-axis for = 1.2. The upper (grid-like) and lower surfaces are for Sm01 
and i^Einst respectively. 

inition also yields the same result agreeing with linear theory [59]) whereas 
Em01 = M — e^/r. This question must be considered important. M0ller's 
energy- momentum complex is not constrained to the use of any particular 
coordinates (unlike the case of the Einstein complex); however, as we have 
shown above, it does not furnish expected result for the Reissner-Nordstrom 
space-time. We agree with Virbhadra's conclusion that the Einstein energy- 
momentum complex is still the best tool for obtaining energy distribution in 
a given space-time. 
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Chapter 8 

Summary and Conclusion 



An important feature of conserved quantities such as the energy, momen- 
tum and angular momentum is that they play a very fundamental role in 
any physical theory as they provide a first integral of equations of motion 
(Nahmad-Achar and Schutz[48]). These help to solve, what would other- 
wise be, intractable problems, for instance, collisions, stability properties 
of physical systems etc. Conservation laws of energy-momentum, together 
with the equivalence principle, played a significant role in guiding Einstein's 
search for generally covariant field equations. Evidently, it is desirable to 
incorporate conserved quantities in general relativity. Energy-momentum is 
an important conserved quantity whose definition has been a focus of many 
investigations. Unfortunately, there is still no generally agreed definition of 
energy and momentum in general relativity. 

Einstein's formulation of energy- momentum conservation laws in the form 
of a divergence to include contribution from gravitational field involved the 
introduction of a pseudotensor quantity tj\ Owing to the fact that t/ is not 
a true tensor (although covariant under linear transformations), Levi-Civita, 
Schrodinger, and Bauer expressed some doubts at the vahdity of Einstein's 
energy-momentum conservation laws. Although, Einstein defended the use 
of a pseudotensor quantity to represent gravitational field and showed that 
his energy-momentum pseudocomplex provides satisfactory expressions for 
the total energy and momentum of closed systems, the problems associated 
with Einstein's energy-momentum complex, used for calculating the energy 
and momentum distribution in a general relativistic system, was followed by 
many definitions, some of which are coordinate dependent and others are 
not. The physical meaning of these was questioned, and the large number 
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of the definitions of energy-momentum complexes only fuelled scepticism 
that different energy-momentum complexes could give unacceptable different 
energy distribution for a given space-time. The problems associated with 
energy-momentum complexes resulted in some researchers even doubting the 
concept of energy-momentum localization. 

Misner et al [45] argued that to look for a local energy-momentum is look- 
ing for the right answer to the wrong question. They further argued that en- 
ergy is only localizable for spherical systems. Cooperstock and Sarracino [22] 
countered this point of view, arguing that if energy is localizable in spherical 
systems then it is localizable in any space-times. Bondi[6] noted that a nonlo- 
cahzable form of energy is not admissible in general relativity. The viewpoints 
of Misner et al discouraged further study of energy localization and on the 
other hand an alternative concept of energy, the so-called quasi-local energy, 
was developed. To date, a large number of definitions of quasi-local mass have 
been proposed. The uses of quasi-local masses to obtain energy in a curved 
space-time are not limited to a particular coordinates system whereas many 
energy-momentum complexes are restricted to the use of "Cartesian coordi- 
nates." Penrose[53] emphasized that quasi-local masses are conceptually very 
important. Nevertheless, the present quasi-local mass definitions still have 
inadequacies. For instance, Bergqvist[5] considered quasi-local mass defini- 
tions of Komar, Hawking, Penrose, Ludvigsen-Vickers, Bergqvist-Ludvigsen, 
Kulkarni-Chellathurai-Dadhich, and Dougan-Mason and concluded that no 
two of these definitions give agreed results for the Reissner-Nordstr0m and 
Kerr space-times. The shortcomings of the seminal quasi-local mass defini- 
tion of Penrose in handling the Kerr metric are discussed in Bernstein and 
Tod[4], and in Virbhadra[74] . On the contrary, the remarkable work of Virb- 
hadra, and some others, and recent results of Chang, Nester and Chen have 
revived the interest in various energy-momentum complexes. 

Virbhadra, and co-workers considered many space-times and have shown 
that several energy-momentum complexes give the same and acceptable re- 
sults for a given space-time. Aguirregabiria et al. [1] proved that several 
energy-momentum complexes "coincide" for any Kcrr-Schild class metric. 
Virbhadra [74] showed that for a general non-static spherically symmetric 
metric of the Kerr-Schild class, the energy-momentum complexes of Einstein, 
Landau and Lifshitz, Weinberg and Papapctrou furnish the same result as 
Tod obtained using the Penrose quasi-local mass definition. These are fasci- 
nating results. Recently, Chang, Nester and Chen [14] demonstrated that by 
associating each of the energy-momentum complexes of Einstein, Landau and 
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Lifshitz, M0ller, Papapetrou, and Weinberg with a legitimate Hamiltonian 
boundary term, then each of these complexes may be said to be quasi-local. 
Quasi-local energy-momentum are obtainable from a Hamiltonian. Hence 
energy-momentum complexes are useful expressions for computing energy 
distributions. 

Virbhadra and Parikh [67] calculated, using the energy-momentum com- 
plex of Einstein, the energy distribution for a spherically symmetrically 
charged black hole in low-energy string theory and found that the energy 
is confined to the interior of the holes. Using Einstein's energy-momentum 
complex, Chamorro and Virbhadra [12] studied the energy distribution as- 
sociated with static spherically symmetric charged dilaton black holes for an 
arbitrary value of the coupling parameter which controls the strength of the 
dilaton to the Maxwell field, and got an acceptable result. We [77, 78], (for a 
discussion see chapter 3) computed energy distributions in these space-times 
using the Tolman form of the Einstein's complex and confirmed both of the 
above computations. In the case of static spherically symmetric charged dila- 
ton black holes the energy distribution depends on the value of the coupling 
parameter while the total energy does not depend on this parameter. 

Earlier investigations [12, 13, 21, 63, 64, 65, 66, 67, 68, 72] with many 
asymptotically fiat space-times indicated that several energy-momentum com- 
plexes give the same and acceptable result for a given space-time. Rosen and 
Virbhadra [56, 71] showed, using the Einstein- Rosen space-time, that even 
for an asymptotically nonfiat space-time many energy-momentum complexes 
could give the same and persuading results. We[79] computed the energy 
distribution in the Ernst space-time, using the Einstein energy-momentum 
complex and got encouraging results. This prompted us to investigate energy 
distribution in Melvin's magnetic universe (which is a special case of Ernst 
space-time) using several different energy- momentum complexes. We [81] 
found that the energy-momentum complexes of Einstein, Landau and Lif- 
shitz, and Papapetrou give the same and acceptable energy distribution in 
Melvin's magnetic universe(for a discussion see chapter 4) . These results 
uphold the importance of energy- momentum complexes. 

The work of Rosen[55] and Cooperstock [18] on the energy of the universe 
was followed by the investigations of Johri et al [37], and Banerjee and Sen [3]. 
These researchers studied the energy content of the universe using different 
models. Using the Einstein energy-momentum complex, Rosen [55] studied 
the total energy of a closed homogeneous isotropic universe described by 
the Friedmann-Robertson- Walker (FRW) metric and found that to be zero. 
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Cooperstock [18] concluded, after making use of calculations involving killing 
vectors, that for a conformal Friedmann-Robertson- Walker metric the total 
energy density is equal to zero. Johri et al. [37] showed, using the Landau 
and Lifshitz definition, that the total energy of an FRW spatially closed 
universe is zero at all times irrespective of equations of state of the cosmic 
fluid. They also showed that the total energy enclosed within any finite 
volume of the spatially flat FRW universe is zero at all times. Using the 
energy-momentum complex of Einstein, Banerjee and Sen [3] showed that 
the energy and momentum density components vanish in the Bianchi type 
I space-times. We [80] extended the investigations of Banerjee and Sen with 
the energy-momentum complexes of Landau and Lifshitz, Papapetrou, and 
Weinberg to check whether these complexes could give different results for 
the Bianchi type I space-times. We got the same results as obtained by 
Banerjee and Sen [3] (for a discussion see chapter 5). 

The Kerr-Newman (KN) solution is the most general black hole solution 
to the Einstein-Maxwell equations. Cohen and de Felice [16] investigated 
energy distribution in this space-time using Komar's prescription. This was 
followed by the investigations of Virbhadra [63, 64] and Cooperstock and 
Richardson [21] who showed (up to the third order and seventh order, respec- 
tively, of rotation parameter) that the energy-momentum complexes of Ein- 
stein and Landau-Lifshitz give the same and reasonable energy distribution 
in KN space-time. Aguirregabiria et al. [1] performed exact computations for 
the KN energy distribution, in the prescriptions of Einstein, Landau-Lifshitz, 
Papapetrou, and Weinberg (ELLPW) in Kerr-Schild Cartesian coordinates. 
They showed that the ELLPW complexes again gave the same energy dis- 
tribution, but the second term of their result differs by a factor of two from 
that obtained by Cohen and dc Felice. Seeing the results that the ELLPW all 
give the same energy distribution, and Lcssner's[41] conclusion that M0ller's 
energy-momentum expression is a powerful representation of energy and mo- 
mentum we found it tempting to obtain energy using this prescription. We 
first note that Florides [26] showed that for all static or quasi-static space- 
times, the M0ller's energy formula is equivalent to the Tolman's energy for- 
mula (Eq. 2.27). We[82] evaluated the energy distribution for the KN space- 
time in M0ller's prescription. We found that the energy distribution in KN 
space-time computed using M0ller energy-momentum complex agrees with 
Komar mass obtained by Cohen and de Fehce[16]. We also found that our 
results support the Cooperstock hypothesis (for a discussion see chapter 6). 

As already discussed, it has been shown with examples of many space- 
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times that several energy-momentum complexes give the same and acceptable 
energy distribution for a given space-time. Recently Virbhadra[74] investi- 
gated whether or not several energy-momentum complexes give the same 
result for the most general nonstatic spherically symmetric metric, and con- 
trary to previous results, he found that the prescriptions of Einstein, Landau 
and Lifshitz, Papapetrou, and Weinberg all give different results. Based on 
some analysis of the results known with many prescriptions for energy dis- 
tribution (including some well-known quasi-local mass definitions) in a given 
space-time Virbhadra[74] concluded that the formulation by Einstein is still 
the best one. In order to test the validity of Virbhadra's conclusion and 
to further investigate the Cooperstock hypothesis we[83] used the M0ller's 
energy-momentum complex to study the energy distribution in the most 
general nonstatic spherically symmetric space-time considered by Virbhadra 
[74]. The M0ller energy distribution differs in general from that obtained 
using the Einstein energy-momentum complex. Both prescriptions agree for 
the Schwarzschild, Vaidya and Janis-Newman-Winicour space-times, but dis- 
agree for the Rcissner-Nordstrom space-time. For the Reissner-Nordstrom 
space-time both the Einstein prescription and the Penrose quasi-local mass 
definition yields the same result which agrees with linear theory. This con- 
firms Virbhadra's conclusion that Einstein's prescription is still the best tool 
for finding energy distribution in a given space-time. 

The main weaknesses of energy-momentum complexes is that most of 
these restrict one to make calculations in "Cartesian coordinates", and the 
large number of these energy-momentum complexes makes it difficult to de- 
cide as to which one to use to compute energy-momentum distribution - given 
the suspicion these would give different energy-momentum distributions for 
a given space-time. The alternative concept of quasi-local mass is more at- 
tractive because these are not restricted to the use of any special coordinate 
system. There is a large number of definitions of quasi-local masses. It has 
been shown[5] that for a given space-time many quasi-local mass definitions 
do not give agreed results. On the other hand previous results [1] and our 
results [79, 80, 81] show that for many space-times several energy-momentum 
complexes give the same and acceptable energy-momentum for a given space- 
time. The important paper of Chang, et a/ [14] dispels doubts expressed about 
the physical meaning of these energy-momentum complexes. Our results 
[82, 83] support Virbhadra's conclusion that Einstein's energy-momentum 
complex is still the best available method for computing energy-momentum 
in a given space-time. These results also support the Cooperstock hypothesis 
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that the energy and momentum in a curved space-time are confined to the 
regions of non- vanishing energy-momentum tensor T- *^ of the matter and all 
non-gravitational fields. 
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